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Introduction. Transverse vibrations of a bimorph consisting of two piezomagnetoelectric layers and located in the 
alternating magnetic field are investigated. Piezomagnetoelectric layers are multilayer composites with alternating 
piezoelectric and piezomagnetic layers. The mechanical and physical properties of such a composite are given by 
known effective constants. 

Materials and Methods. The applied theory of multilayer plate vibrations takes into account the nonlinear distribution 
of electric and magnetic potential in piezoactive layers in the longitudinal and transverse directions. On the basis of this 
theory, the stress-strain state, the dependences of deflection, electric and magnetic potentials on the volume ratio of the 
composition of the hinged bimorph, are investigated. The electric potential 1s assumed to be zero at all electrodes, while 
the magnetic potential is zero at the inner boundary and unknown at the outer boundaries. Therefore, the distribution of 
electric and magnetic potentials in the middle of the layer are unknown functions. In the case of the magnetic potential, 
the distribution at the outer boundary is also unknown. In the problem, the Kirchhoff hypotheses for mechanical 
characteristics were accepted. The use of the variational principle and the quadratic dependence of the electric and 
magnetic potentials on the thickness of piezoactive layers made it possible to obtain a system of differential equations 
and boundary conditions. 

Results. When the volume ratio of the composition of piezoactive bimorph materials changes, the electric potential in 
the middle of the layer changes nonlinearly. The magnetic potential in the middle of the layer and at the outer boundary 
increases almost linearly with an increase in the volume percentage of Ba7iO3. The dependence of the deflection in the 
middle of the layer is determined. 

Discussion and Conclusions. An applied theory for calculating transverse vibrations of a bimorph with two 
piezomagnetoelectric layers is constructed. The dependence of the characteristics of the stress-strain state, electric and 


magnetic fields on the volume fractions of piezomagnetic and piezoelectric materials, is investigated. 
Keywords: piezoelectrics, piezomagnetics, composite, bimorph, magnetoelectroelasticity, bending vibrations. 
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Introduction. In the production of sensor and measuring systems, small household appliances, cell phones, 





and wireless sensor systems, powerful energy sources are not needed to monitor and diagnose the technical condition of 
objects. In this case, the prerequisites are the mobility and volatility of the above devices. 

Piezoelectric materials directly convert electrical energy into mechanical energy and vice versa. This property 
allows them to be widely used in science and technology. These materials are used in ultrasonic emitters of elastic and 
acoustic waves, receivers of such waves, devices for suppressing vibrations of machine elements and structures, etc. 
Recently, another field of application of piezoelectrics has been rapidly developing — energy collection and storage 
devices. In this case, piezoelectric materials are part of piezoelectric energy generators (PEG). PEG are placed on 
elements of machines or structures that vibrate intensively, which are in the zone of elastic wave propagation or are 
exposed to variable pressure. The main types of these devices have a bimorph or stack multilayer structure and 
experience bending or longitudinal deformations, respectively. Low-power sources of electric current are created on the 
basis of PEG. They include autonomous power sources (e.g., for damage monitoring devices in hard-to-reach places of 
pipeline structures, etc.). An overview of such devices is available in [1—2]. One of the ways to design effective PEG is 
the use of piezoactive composites of various types of connectivity and heterogeneous materials based on piezoceramics, 
including porous one. 

PEG, in whose design there are additional electromagnetic elements or permanent magnets, can fix or use the 
energy of an alternating magnetic field. One of the ways to solve this problem is the use of piezomagnetic materials in 
combination with piezoelectric ones. In this case, the alternating magnetic field causes the deformation of the 
piezomagnetic and the coupled piezoelectric, as a result, the latter generates electrical energy. There is a class of 
materials with ferromagnetic properties. Piezomagnetism is a phenomenon observed in some antiferromagnetic and 
ferromagnetic crystals. It is characterized by a linear relationship between the magnetic polarization of the system and 
mechanical deformation. In a piezomagnetic material, a spontaneous magnetic moment can be induced through 
applying mechanical stress, or deformation by applying a magnetic field. In studies of piezomagnetic materials [3-5], 
CoFe,O41s very often considered. In [6-8], a composite based on CoFe,O, and BaTiO; with piezoelectric and 
piezomagnetic properties is studied simultaneously. 

Solutions to the problems of electroelasticity and magnetoelasticity are given in [9-11]. In [12], applied 
theories of vibrations of multilayer piezoelectric plates were developed considering the specifics of the distribution of 
electric potential over the structure thickness. 

The problems on steady-state oscillations of an electro-magneto-elastic layer and a half-space under the action 
of harmonic loads are presented in [13, 14]. The prestressing is taken into account, as well as various electrical and 
magnetic conditions at the boundaries. The effect of these factors on the dispersion properties is investigated. 

Earlier [15, 16], an applied theory was developed that considers the inhomogeneous distribution of the electric 
potential in the longitudinal direction, and the quadratic dependence on thickness. In the same papers, the stress-strain 
and electrical state of a hinged and cantilevered bimorph is investigated. In both cases, the applied theory showed good 


convergence with the finite element modeling results. The authors also developed an applied theory of bimorph 
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vibrations [17] consisting of an electroelastic and magnetoelastic layer. This approach is in good agreement with the 
finite element analysis results. 

In this paper, the vibrations of the device are considered in the framework of a plane deformation. Based on the 
variational principle, an applied theory of bending vibrations of a two-layer piezoelectric bimorph is constructed. For 
steady-state vibrations, boundary conditions and a system of differential equations are obtained for four unknown 
functions (deflection, electric potential in the middle of the layer, magnetic potential in the middle of the layer and at 
the outer boundary), depending on the length of the bimorph. The influence of different percentage volume ratios of the 
bimorph composition on deflection, electric and magnetic potentials in certain positions, is investigated. The study 
results provide selecting the makeup of a composite piezomagnetoelectric material to achieve the most efficient 
operation of the device. 

Materials and Methods. A plate consisting of two identical piezoelectric layers is considered. It performs 
steady-state transverse vibrations within a plane deformation. Each layer is a 2—2 connectivity composite consisting of 


alternating piezoelectric and piezomagnetic layers (Fig. 1). 








Piezoelectric layer =a Piezomagnetic layer 


Fig. 1. 2-2 composite structure 
Effective properties of such a composite were found in [8]. Large surfaces of the layers are electrodated, and 
the layers themselves are polarized in thickness. The bimorph is hinged at the edges, all surfaces are free from 
mechanical stresses. The upper and lower boundaries of the plate are affected by magnetic flux B,, while at the 
boundary between the layers, the magnetic potential 1s considered to be zero. The electrical potential is zero at all 
electrodes. The side surfaces are considered to be insulated from magnetic and electric fields. 
The equations for describing the vibrations of a composite with effective properties, connectivity of 
mechanical, electric and magnetic fields, have the form [18]: 
V-o+pf =pu, V-D=0,, V-B=0, 
G6=cit-e -B=ir -H, 
D=e:6+k-E+a-A, 
B=h:e+a' -E+u-H, 


e=—(Vu+(Vu)'), Hao, Fs. (1) 
Here, O and € — mechanical stress and strain tensors, D and E — vectors of electric induction and electric field 
strength, B and H — vectors of magnetic induction and magnetic field strength, p — material density, c — elastic 
moduli tensor, e — piezoelectric moduli tensor, / — piezomagnetic moduli tensor, « — dielectric permittivity 
tensor, o& — magnetoelectric moduli tensor, 4 — magnetic permeability tensor, / — mass force density vector, 6, 
— electric charge volume density, u — displacement vector, @ and € — electrical and magnetic potentials. 


The boundary conditions are determined for the mechanical, electric and magnetic fields, respectively. 


For the first case, we note the absence of mechanical stresses at the bimorph boundary: 


o,,-7,| 


The biomorph is hinged at the ends (Fig. 2): 


UU, f7 13: 
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Fig. 2. Geometry and boundary conditions of bimorph with composite piezomagnetoelectric layers 


Next, we formulate the electrical boundary conditions. Electrical potential on the internal and external 
electrode, respectively: 


| =f = Vi. Ol, i _ V, 7 


We indicate the magnetic boundary conditions. Magnetic potential at the inner boundary: 
é|  =M,. 


x3 =0 


Magnetic flux 8B, affects the upper and lower boundaries of the plate: 
=B 


I. =tH 0° 
We use the variational equation for steady-state vibrations [10]. It generalizes Hamilton's principle in the 
electroelasticity theory taking into account magnetic components. For the case of plane deformation in the absence of 


surface loads and in the presence of magnetic flux: 


| | SHdS — pw’ | | u,du.dS + | | (p8u, +,59 + B,dE)dS =0, (2) 
S S Ss 


where 5H = Ogee D,6E;- B,5H;. 
To construct an applied theory of vibrations, we will accept Kirchhoff's hypotheses. In accordance with them, 


the distribution of displacements along the thickness has the form: 


Uy (X1,X3) =-X3Wi1, Ug (X1,X3) = W(X) (3) 

In particular, the single normal hypothesis is accepted for a mechanical field. Next, we consider a problem in 

which the value of the electric potential on the electrodes can be zero, so its distribution is not described by a linear 
function. Taking into account the possible inhomogeneity in the length of the element associated with the influence of 


boundary conditions at the ends of the bimorph, its thickness distribution is assumed to be quadratic: 
Ma 4X; Ne 2X 
x,,X,;) =V,(x,) =| —-1]4+V, (x +V. <=} +41]. 4 
(x5) = so) [ 7 {- ls (x) =| 2 +1) (4) 


Here, X; = x3-h/2. Functions Vo, V, and V, are responsible for the value of the electric potential at the inner electrode, 
in the middle of the layer, and at the outer electrode, respectively. To satisfy the conditions of the problem, we take 
these functions in the following form (see Fig. 2): 
Vi(~%) =Y =const, V(x) =P), V(x) =V, = const. 
Here, function ®(x,) is unknown. 
We represent a quadratic distribution of the magnetic potential over the thickness of each layer. The 


distribution along the length is heterogeneous, at the inner boundary of the layers, its value is assumed to be zero: 


(x,,%,) = Mone [ ae eam cxp[- + Jeaseay3(22-1). (5) 
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Here, X3 = X3-—h/2. Functions My, M, and M, are responsible for the value of the magnetic potential at the inner 
boundary, in the middle of the layer, and at the outer boundary, respectively, and are taken as follows (Fig. 2): 
M,(~%) =M, =const, M,(x%,)==,(%), M,(~%) =2,(%). 
Here, functions £,(x,) and &3(x,) are unknown. 
We substitute relations (3)-(5) into equation (2) and integrate it by the bimorph thickness, and then we 
equate the coefficients with independent variations 5w, 6®, d=, and d#, to zero. Thus, we obtain a system of four 
differential equations (6) from four unknown functions depending on x, (then, we omit the subscript), and five 


boundary conditions (7). 


















































a aa ae aE 2 ee (x)- 
pean Phan E(x) + _— ie w(x) = 2B, -E# M, =0, 
= ca (x) + “2 - cme: 9 sh “ ,(x)-2p,H - 
3 __ (6) 
Bet yea pu(x)H + _ (x) =0 
ra a 
FOO BROS gnome 
fully, Sully, Bull gy), Mulls (yy Stull s (4 2M A iy) all og, =o, 
a D(x) ill ds 5, (x) + all ds (yy 2H Er dy), POM be), (7) 


* ‘ ‘ ‘ = o) ed ~ 
Here, the following designations were introduced: ¢, =G,—-¢3/C33, €31 = €31-C13€33/C33, h31 = 


ha,- C13h33/C33, O33 —— X%33- C33h33/C33, E33 = Con Oss" (Cae ° They occurred after satisfying condition 033 — 0 and 
exclusion of €33. 
Research Results. The results of the bimorph calculation according to the proposed theory are compared to the 


finite element calculation in the low-frequency region for the volume ratio of the piezoelectric and piezomagnetic 
components 80 % BaTiO; and 20 % CoFe 0,4. The comparison has shown that the error in finding the characteristics of 
the mechanical and magnetic fields is less than 1%. When determining the electric field in the middle part of the plate, 
the difference was about 5%. Describing the situation in the vicinity of the support points, it should be noted that the 
size of the neighborhood along the longitudinal coordinate 1s approximately equal to the bimorph thickness. Here, when 


determining the electric field, a difference of 20% is recorded. 
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The first step in studying the vibrations of a two-layer piezomagnetoelectric bimorph with a change in the 
volume ratio of BaTiO; and CoFe,O, in the composite is to determine its effective properties. Tables 1 and 2 present 
these properties found from the results of study [8]. 


Table 1 
Material constants (elastic modules, dielectric and magnetic permittivity) 
for different BaTiO; volume fraction 


Volume Elastic modules Dielectric permittivity | Magnetic permittivity 


fiction [GPa 10m TON STC 
coat ee eee ee 


0 | 2860 | 1730 | 1700 | 2695 | 45.30 | 0080 | 0.093 | 5.900 | 1570 
10} 2709 | 160.4 | 1549 [260.0 [45.07 | 1469 [0073 | 5318 | 0.632 


6 | 2067 | tos.t | 102.1 | 2156 | 4392 | 7490 | 0.220 | 2.390 | 0.155 
80 | 1860 | 923 | 859 | 2013 | 4346 | 9.488 | 048i 1200 | 0.120 
99 | 1766 | 854 | 789 | 1989 | 4323 | 10353 | 0.857 | 0635 | 0.110 


Table 2 





Material constants (piezoelectric, plezomagnetic and magnetoelectric modules) 
for different BaTiO; volume fraction 


Volume 
A A, SSS 
BaTiO; (%) eS SS eee ee 
of 
ea [ores ee es || 


x0 2a 


- | 0038 | 0220 [696 | 297 [228 | 20 | 34s | 2a 

80 | 0060 [057 [928 [137 91110 
[90 | 0263 | 1e7_| roa [46 [a6 ‘| 55 ([ 2 
[100 [4.400 | is600 [eo fo fo jo [0 | 0 





The bimorph vibrations were excited by a magnetic flux applied to the upper and lower faces (Fig. 2), which 
varied according to the harmonic law with an amplitude of B, =5x10~° Wb and a frequency of 10 kHz. 

Figure 3 shows the deflection in the middle of the layer depending on volume fraction of BaTiO. It can be 
seen from the graph that the deflection in the position having coordinates x, =L/2, x, =H/2, 1s zero if the bimorph 
consists only of a piezoelectric BaTiO;. The deflection of the bimorph reaches the greatest value if it contains only 


piezoelectric magnet CoFe,O,4. The deflection almost linearly depends on the volume ratio of the components of 


piezoactive materials. 
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Fig. 3. Deflection w(x, ) in the middle of the layer for different BaTiO; volume fraction 


Based on the data in Figure 4, it can be concluded that the electric potential in the middle of the layer varies 
nonlinearly with a change in the volume ratio of the composition of bimorph piezoactive materials. If the bimorph 
consists only of BaTiO; or CoFe,O,, then the electric potential at the point (L/2,H/2) is zero and reaches the highest 


value at 35 % BaTiO; in the bimorph. 
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Fig. 4. Electric potential ‘Y(x,) for different BaTiO; volume fraction 


The analysis of Figures 5 and 6 allows us to conclude that the magnetic potential in the middle of the layer =, (L/ 2) 


and at the outer boundary =,(L/2) increases almost linearly with an increase in the volume of BaTiO; in the bimorph. 
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Fig. 5. Magnetic potential &,(L/2) for different BaTiO; volume fraction 
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Fig. 6. Magnetic potential &,(L/2) for different BaTiO; volume fraction 


Discussion and Conclusions. An applied theory is proposed for calculating transverse vibrations of a bimorph 
made of two layers of a composite based on CoFe,O, and BaTiO; with both piezoelectric and piezomagnetic properties, 
in an alternating magnetic field. Such a design can serve as a model of a piezoelectric generator of a device for 
collecting and storing energy under the action of an external magnetic field. In the low-frequency region (below the 
natural frequency of the first bending mode), calculations of the stress-strain state of the bimorph, the distribution of 
electric and magnetic fields, are carried out. The dependence of the deflection, electric and magnetic potentials on the 
volume ratio of the bimorph composition is investigated. In further work, it 1s assumed to determine the output potential 
and power of an electric current excited by an alternating magnetic field. The purpose of these surveys will be to collect 
electrical energy. 
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Introduction. The paper considers an axisymmetric problem of elasticity theory for a radially inhomogeneous 
transversally isotopic nonclosed sphere containing none of the 0 and z poles. It is believed that the elastic moduli are 
linear functions of the radius of the sphere. It is assumed that the side surface of the sphere is fixed, and arbitrary 
stresses are given on the conic sections, leaving the sphere in equilibrium. The work objective is an asymptotic analysis 
of the problem of elasticity theory for a radially inhomogeneous transversally isotropic sphere of small thickness, and a 
study of a three-dimensional stress-strain state based on this analysis. 

Materials and Methods. The three-dimensional stress-strain state is investigated on the basis of the equations of 
elasticity theory by the method of homogeneous solutions and asymptotic analysis. 

Research Results. After the homogeneous boundary conditions set on the side surfaces of the sphere are met, a 
characteristic equation is obtained, and its roots are classified with respect to a small parameter characterizing the 
thickness of the sphere. The corresponding asymptotic solutions depending on the roots of the characteristic equation 
are constructed. It is shown that the solutions corresponding to a countable set of roots have the character of a boundary 
layer localized in conic slices. The branching of the roots generates new solutions that are characteristic only for a 
transversally isotropic radially inhomogeneous sphere. A weakly damping boundary layer solution appears, which can 
penetrate deep away from the conical sections and change the picture of the stress-strain state. 

Discussion and Conclusions. Based on the solutions constructed, it is possible to determine the applicability areas of 
existing applied theories and propose a new more refined applied theory for a radially inhomogeneous transversally 
isotropic spherical shell. 


Keywords: equilibrium equations, Legendre equations, radially inhomogeneous sphere, characteristic equation, 
boundary layer solutions, variational principle, applied theory, reduction method. 
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Introduction. One of the properties of materials affecting the stress-strain state of elastic bodies is their 
heterogeneity. Investigating the stress-strain state of inhomogeneous bodies based on three-dimensional equations of 





elasticity theory is associated with significant mathematical difficulties. 

A number of studies are devoted to the investigation of three-dimensional problems of elasticity theory for the 
sphere. 

In [1], based on the equations of elasticity theory for the sphere, a general solution satisfying the boundary 
conditions on the contour in the sense of Saint-Venant was obtained, the stress-strain state of the sphere was analyzed. 
In [2], based on the equations of elasticity theory for a thick isotropic sphere, homogeneous solutions depending on the 
roots of the transcendental equation are constructed. In [3], on the basis of solving three-dimensional problems of 
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elasticity theory for a sphere of small thickness, the accuracy of existing applied theories is studied, and a method for 
constructing refined applied theories is given. A three-dimensional asymptotic theory of a transversally isotropic 
spherical shell of small thickness is presented in [4]. An analysis of the three-dimensional stress-strain state of a three- 
layer sphere with a soft filler is described in [5]. In [6], the torsion problem is studied by the method of homogeneous 
solutions for a radially inhomogeneous transversally isotropic sphere of small thickness, when the elastic characteristics 
are changed by linear, quadratic and inversely quadratic laws along the radius. In [7], the torsion problem for a radially 
layered sphere with an arbitrary number of alternating hard and soft layers is studied. The existence of weakly damping 
boundary-layer solutions and a possible violation of the Saint-Venant principle in its classical formulation are shown. 
An applied theory of torsion of a radially layered sphere is constructed, which adequately takes into account the 
emerging features. In [8], the problem of elasticity theory for a radially inhomogeneous hollow ball is investigated using 
the finite element method and spline collocation. In [8], using the finite element method and spline collocation, the 
problem of elasticity theory for a radially inhomogeneous hollow ball is studied. The results obtained by finite element 
methods and spline collocation are compared. The axisymmetric problem of the theory of elasticity for a radially 
inhomogeneous transversally isotropic sphere of small thickness is studied by the method of asymptotic integration of 
the equations of the theory of elasticity in [9]. Inhomogeneous and homogeneous solutions are constructed. The nature 
of the stress-strain state is established. In [10], an axisymmetric problem of elasticity theory for a sphere of small 
thickness with variable elasticity moduli is considered by the method of homogeneous solutions. Asymptotic formulas 
for displacements and stresses are obtained, which provides calculating the three-dimensional stress-strain state of a 
radially inhomogeneous sphere. 

Materials and Methods. Deformation is considered within the framework of the linear theory of elasticity of a 
nonclosed sphere, whose material is transversely isotropic and inhomogeneous along the radial coordinate. The 
thickness of the hollow sphere is assumed to be small compared to the radius and size along the arc coordinate. 
Boundary conditions are considered that provide solving the problem in an axisymmetric formulation. We assume that 
the sphere does not contain any of the poles QO and z. In the spherical coordinate system, the area occupied by the 
sphere will be denoted by [={re[n;7], 9¢[6; 4], d€[0; 2z]}. 

The linear dependence of the elastic — of the material along the radius is considered: 

A, =a)'r, A, =ayr, A, =ayr, A), = ay, Ay, = Ayr, i) 
where ai)’, a\;”, as) ,asy’, ay’ — some constants. 

The system of equilibrium equations in the absence of body forces in a spherical coordinate system r,0,@ has 
the form [11]: 

Oo, n 100,, 7 26, Oy Og +C,,Cled 





'r — 0, > 
or r 0o@ r ) 
00, 1 OO pg : 30,9 + (Ogg — Tyg )ctg@ 0 3) 

or r oO r 


where, 0,,.,0,9,O 5 


Dv. =v.(r,9@), D0, =v, (7,9) as follows [4]: 


O,, — Stress tensor components, which are expressed in terms of displacement vector components 















Ov. A, Ov 
= r+) py cte94+2v, +— |, 5 
rr 11 ar y ( 0 S r ve ( ) 
O45 = A, | (A, + A, )v, + A,,Uz dale 5 (6) 
| Ov, 
Ogg = Ar +—| (A, +A )0, +A,,0,¢120 + A, — |, (7) 
Or Ir (ale, 
Ov, OU, 10D 
O., = A,,| —>-—++-— |. 8 
ro ul ar y 7 . ( ) 


Substituting (5)-(8) into (2)-(3), taking into account (1), we obtain the equations of equilibrium in 
displacements. 
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Ou Ou Oru Cru. Ou 
xe" | —24+u,cted |+e(b +69 )| —2 cte@ + |+ 6*bO | —“+— cite |=0, 9 
Se gS (d:; 12 le § a |S 00 § (9) 
: Cru O 
(| SMe 4 26 Me 364, +e |+ (300 +5 +6 Je? “24 
Ap Ap p00 a0 
21(0) | Og OU » (0) Ou, 2 (1.0) 2.2 (0) 
+2"by) | cig +— > | + eb ao (d') ctg*O + bY Ju, =0. (10) 


eA i 


1 ; ; 1 int 
Here, p= Lin{ =) — new dimensionless variable; ¢ = sn( 2 | — a small parameter characterizing the thickness of 
0 


(0) 
a, 


B) D é : - 
the sphere; 7, = Nic b. pel-i,1), v= — Uy = ? bi ) = a dimensionless quantities; G, — some parameter 
0 0 0 


having the dimension of the elasticity modulus. 
We assume that the lateral part of the sphere boundary is fixed, 1.e., 


=, (11) 
Ug| 1, =0- (12) 
We assume that at the ends of the sphere (on the conical sections), the stresses are specified 
Cvleo, = Ss)» — Fpaly.g = Sx(P)- (13) 
Here, f,,(?), f5,(P) (s =4.2)— sufficiently smooth functions that satisfy the equilibrium conditions. 


Solutions (9), (10) are sought in the form [3, 4]: 
u,(p,8)=a(p)m(A); u,(p,A) =d(p)m'(8), (14) 


Where function m(@) satisfies the Legendre equation: 
m"(O)-etg0-m'(O)+{ 2 5 me =. (15) 
After substituting (14) into (9), (10), (11), (12), taking into account (15), we obtain: 


" f 1 
Hi a"(p) + 2eti%a"(p) +6 (an -26 2080) [22-15 |a(p)- 


(16) 
I 
-[2-j e[(62 +A )a'(o)- 0(0 +02 +08) - 269) acy] =0. 
bi! (d"(p) +2ed"(p))-e" (2 1 ja +(B -B© +360 Jaw + 
(17) 
2° (30.0 +80 +00 Ja(p)+e (BS +8 Jap) =0. 
a(p) =0, npu p=. (18) 
d(p)=0, upu p=. (19) 
The solution to system (16), (17) has the form: 
a(p) =e” pe" C, + pe ’C, + p,e”"C, + pe ”C, | , (20) 
d(p)=e" | 4c" C ge" Che €.7¢.6 "CC | ; (21) 


where C,, (n =1,4) — arbitrary constants, 
p=), s. -(# + Joe +(BS) — BS — 4p") ) ; 


_ (0) (0) (0) (0) (0) (0) \. 
ty =—(bf +B )s, —(258 +b, +b;; — bi i; 
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q, = (a{) +b )s, —(25 + bY +b) —bY ); s, — equation roots 
1 
(0) 7,(9) 4 2 (0) \2 (0) 7,(0) (0) 7,(0) (0) 7,(0) (0) 7,(0) (0) 7,(0) 
by Oya 8 (: 4 ((d); y+ 2b By) — yy by, )+5{ by — Sb byy + 4b,y by’ — 


1 
-26H! — HH — 2H! ] 2 4 (2 : + os ~ 4b +26 +260 +a) (22) 


l l 
«(2 + Joe +b —B +400 | -(# 5 |(2n8 +B +5 9) =0. 
The system of linear algebraic equations with respect to C,,C,,C,,C,, is obtained through satisfying the 


homogeneous boundary conditions (18), (19). The zero equality of the determinant of this system is a condition for the 
existence of nonzero solutions and leads to a characteristic equation with respect to spectral parameter 
rN C2 é) = (PO = DG, \(t,P, =i), ) sh? (e(s, +S, )) + 03) 
+( pt, — Px, Va — Pyt, ) sh? (e(s, =5, )) =(), 
Equation (23) has countable set of roots Z,. The general solution to the problem is obtained by summing over 


the roots of equation (23) 


u, =>M,a,(p)m, (9), (24) 
u, =>M,4,(pym (0), (25) 


where 
nip=e" | D, (ea, =e" ) + Py (e°""a, —e "a, )| ; 


d.(p)=e” ae O=Oe°" @,+1,e ~ Oj =q,e°"" 0, |, 
eS} =8( 5; 4295 o( 25555 
O, = Pr (4, —, Je —t, (P.4> — Po Je +d, (pt, —~ Pr Je . 


ae +255 ) 


Q, = Pol, (q, —t, je“ —t, LCE — Pyt, jer +> (pt, =.) 


—e(2s, +5, ) 


a t; (Pd —~ P2Q Je —~q (Dd — Pot, jer) + P\Q> (q, ft, ye ’ 


(s5 —25,) ES 


QO, =, Loe — PQ Jee —q (pt, — Prt, er + pyt, (q, fo, je *. 


The set of roots of equation (23) at ¢ +0 consists of countable sets of roots 


Ze ou +O(E). (26) 
E 
For 6),, we have: 
1°. At b >0, b?—b, >0: 
(s, —s, ) (di) —b” ss, )sin((s, +85, )6) +(s, +5, y(n? +h s,s, )sin((s, =S, )d) =0, (27) 
where 
s,= Vb +b =),3 S.= Vb — fb = 
b, = (2620) (26940) +2)? —APBY); by = LY" 
2°. At b >0, b —b, <0: 
B| 26a? —(b)) (a? - B*)-b4) |sh( 25a) + 
(28) 
+a | 24)? p? +(b) (a? - 6?) )] sin (258) =0, 
where 


S; = b, +b; -2, =+(a+if); 
‘= b, —\[b; -b, = +(a —-if). 
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3° At bh >0, b =b,: 
(d\)’s° — Bi? )sin(2sd)+2(Bi)s° +B{) )ds =0, (29) 
where s = Jb, ; 
4° At b, <0, b —b, >0: 


Clay Nr —b” ss, )sh ((s, +85) 6) (Sirs, (ay? +b ss, )sh ((s, =5, )6) =0, (30) 


where 
s, = 16 |-Je? -4,; s, = y\0,|+ fo? 2. 
5°. At b <0, b? —b, <0: 
B| 2b\))a? -(b)) (a? - B*) 04) |sin (25er) + 


(31) 
+a BS p? +(B) (a? - B)-b ) |sh( 268) =0, 
where 
s, = lb |-l2? -2, =+(@-i); 5) = y || + lb? -b, =+(a +i). 
6°. At bh <0, b =b,: 
(d\'s° — Bi )sh(2sd)+2(B\)s° +i) )ds =0, (32) 


where s = ,/|b|. 


Equation (27)—(32) has a countable set of solutions. 
Let us present an asymptotic construction of solutions corresponding to different groups of roots of the 
characteristic equation (23). Substituting (26) into (24), (25) and expanding the resulting expressions in powers € , we 


have: 
Is 
a) u, )= VE: (0? +b )| (a2 +b\?'s; )s, cos(6,s, )sin(d,s,2) 
gol 
—(b) + bs? is cos(6,5; )sin(6,s,p) | + O(e)| m,(@), (33) 


0) = E65: {(22 +bhs (a +b s >) cos( (5,5, )cos (6,8, 0) —cos (d,s, )cos(6,s,2) |+O(e )m (9), (34) 


where 6,, are the solutions to equation 


(s, —s,) . —b” ss, wee 85 )6) +(s, +5, (ay? +b ss, )sin((s, = 5, )d) =0 (35) 
6) uw, = Dae (os be +B))| (BY? +B\s; Js, sin(6,5, )cos(d,s,0)— 
—(b,2 +?'s; )s, sin (5,:, cos (5,s,p) |+O(e)}m, (8). (36) 


Uy (p:0) = ES: {(28? + DOs al (sy $ bs ; )[ sin( (5,5, )sin(d,s,e)—sin(6,s, )sin(d,s,0) |+ O( e) m,(@), (37) 


k=l 


where 6,, are the solutions to equation 


(5, —s,)(BY - —b”s s , )sin((s, +s,)5)—(s, +5,)(BQ +4{)s,s, )sin((s, —s,)6) =0. (38) 


7° 
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a) u,(p;0) = SES: ee +b ){[ Bsin (6, Be) sh(5,ap)— cos (6, 8) ch(5,a@p) | 


| 2086!) cos (5,8) sh(d,a)+(B (a? = B’)-bY )sin(6, B)ch(d.@) = 
-[asin(5, 8p) sh(5,ap) + Bcos(6,h2)ch(5,ap) | 


| 2a@Bb\ sin(5,)ch(6,a)—(B? (a? — B")—b{) )cos(d, B)sh(6,a.) | i ave) m,(0), (39) 


u, (p39) = S65! | 20488 sin(6,8e)ch(6,ap) -(B9 (a? — p’)-b) )x 


xcos(5,A2)sh(5,ap) |] 2a@Bb\? cos(5,8)sh(5,2)+(A? (a? —B")—b,? )sin(6,B)ch(5,a) |- 
-| 206? cos (6B, p) sh(d,ap) +(4. (a? - B*)—bi? )sin(6, Bp) ch (6,ap) x 
| 2aBb\” sin(5,8)ch(5,a) — (a (a? - B’)-BY cos (6, B)sh(6,c) + o(e)} m| (0), (40) 
where 6, are the solutions to equation 
B 26 (.a" — (8? (a? - B?)-b )] sh(26a) +a | 25? B +(BY (a? - B?)-bo )] sin(258)=0. (41) 


ee) 


6) u, (p30) = > BO (BY +b ){) Bsin(5,Bp)ch(6,ap)— acos( 6, Bp) sh(5,ap) |: 


k=] 


| 2a8b\) cos (6,8) ch(5,a) +(b)) (a? - B*)—bY )sin(6, f)sh( 5,2) - 
-| asin(5,Bp)ch(6,ap)+ Bcos(d,8) sh(d,ap) || 2aBb\) sin (6,8) sh(d,a)- 
-(B{? (a? ~B°) -bi) Jeos(5,8)ch(5,a) |+O(e)}m, (8), (42) 
u, (39) = aL | 2x88)? sin(d, 6p) sh(d,a@p) -(B) (a? — B’)—bi! )x 
xcos (5,2) ch(5,ap) || 2a@Bb\? cos (5,8) ch(5,2)+(d? (a? — B”)—b.) )sin(6,8)sh(d,@) |- 
-| 2a}? cos (6B, p)ch(5,ap) + (A? (a? — 6") —b.) )sin(6, fp) sh(d,ap) x 


| 2B” sin(5,8)sh(5,a) -(a? (a? - B’)-BY cos (6, B)ch(5,2) | + o(e)} m| (0), (43) 


where 6,, are the solutions to equation 
B| 20).@ - (A? (a? -B°)-o )] sh(26a)—a | 20)? B +(B9 (a? -B°)-1 )] sin (268) = 0. (44) 
3°. 


a) u,(p;0)= De oes + bi? )(cos (6,8) cos(d,sp) + psin(6,s)sin(d,sp))+ 
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(0) 2 (0) 
i (Di =O, 


SO; 


in(asdon(ase}O(0) m (0), (45) 


: BOs? +b) 
ae ES ( i 44 
“o(e:8) 2 © | (Bs? +5 )5,s 


x(psin(6,s)cos(d,sp)—cos(d,s)sin(d,sp))+ O(e)} m,(@), (46) 
where 6,, are the solutions to equation 


(d\)'s° —bi? )sin(2sd)+2(b))s° +i) )ds =0 (47) 


2g {(B10's? B22) 
6) u,(p;0)= > E,” —.—eos(4,5) sin(5,59) —(8's? + bf )x 


k=1 Ks 


x (sin (6,5) sin(6,s)+ ecos(d,s)cos (5,52)) + O(e)} m, (0) (48) 
co ie +b ) 
u,(p;0)= > ES 4+—_—* x 
099) pe © | (00 +B9)6,s 
x| pcos(d,s)sin(d,sp)—sin(d,s)cos(d,sp) |+ O(e)} m,(@), (49) 
where 6,, are the solutions to equation 
(d\)’s° — bi? )sin(2sd)—2(B))s° +i) )ds =0 (50) 


4°. In case b, <0, b> —b, >0, asymptotic formulas for displacements are obtained from (33)—(38) through 
replacing s,,s, with is,,is,. 

5°. In case b, <0, b> —b, <0, asymptotic formulas for displacements are obtained from (39)—-(44) through 
replacing s,,s, with is,,is,. 

6°. In case b, <0, b> =b,, all asymptotic formulas for displacements are obtained from (45)-(50) through 


replacing s withis . 
For roots (26), the main term of the asymptotic solution of equation (15) at e + 0 takes the form [9, 10]: 


l l : 
= Teo -8 _g2, (0-6, ) (+0); near 0=6, 
v Y Ok 





m,(@) = (31) 
= 0x9 &"' \-5;, (0-8,) (14+ O(6)); near 0 = 0, 
ane 0 [52 Ok 2 2 
We represent displacements in the form: 
Uu, (p, 0) = ye (p)m, (0), (52) 
k=l 
Ug (p, 0) = > Ea, (p)m, (0). (33) 
k=] 
We represent stresses o,, and o,, in the form: 
Oy = > E, (of (p)m, (0) + 0? (p)m, (A)ctgd), (54) 
k=l 


Akhmedov N. K., et al. Analysis of the stress-strain state of a radially inhomogeneous transversely isotropic sphere with a fixed side surface 





é,-) Fo, ©), (55) 
k=1 
here, 
lias | 
ai0(p)=2[APai(p)+0(02 +0.)a,(o)-eb.? (22-5 }a,(0)} 


oi?(p) = (BBS), (p) 
(0) 


b 
O>,(P) = ~ lh (p)+e(a,(p)—d,(p))]- 


The character of solutions (33)-(50) depends essentially on the type of roots 6,,. The first boundary-layer 
terms of these solutions correspond to the Saint-Venant edge effect [4]. In the case of imaginary roots 6,,, these 


boundary layers have weak damping. Thus, the stress-strain state 1s far enough away from the ends and significantly 
depends on them. That is, in this case, the transversely isotropic properties of the inhomogeneous material significantly, 
in comparison to the isotropic material of the sphere, change the pattern of the stress-strain state. At the same time, for 


real or complex 6,,, the pattern of the stress-strain state of an inhomogeneous sphere for such materials qualitatively 


coincides, differing in the decay rate of the above-described Saint-Venant boundary layer solutions of an 
inhomogeneous plate. 


From (51), it turns out that when moving away from the conic sections 0=6,(j =1,2), solutions(33)-(50) 
decrease exponentially. 

Since the constructed solutions satisfy the equilibrium equation and boundary conditions on the side surface, 
the Lagrange variational principle takes the following form [4, 11]: 


YI [los = Fij (p)) du, + (Oyo fj (p)) ou, | | > e?dp=0. (56) 


j=l v= J 
Substituting (52)-(55) into (56) and taking OE , aS independent variations, we obtain an infinite system of linear 


algebraic equations 
oo Ga, (57) 
k=] 


Here, 


1 , ‘ 
dn = Joo (rved pf Sm, (0,) 08.) + food (red p» 
= s=l a 


s=] s=l1 


«(mie ym' (0, re6, fon(p)a, ire*dp| Sm (9,)m, @))} 


2 


= > (0,)| fi(p)d,(p)e”dp+m,(0,)| f,,(P)a, inerdp| 


s=l 

System (57) is always solvable under physically meaningful conditions imposed on the right side (57). The 
solvability and convergence of the reduction method for (57) is proved in [12]. 

Using the smallness of parameter &€ , it is possible to construct asymptotic solutions of system (57). 

Research Results. The structure of the stress-strain state of a radially inhomogeneous transversally isotropic 
sphere of small thickness is analyzed under kinematic conditions on the side surface. It is shown that, in the case of 
fixing the side surface, the character of the solution is determined by the boundary layers. It is found that the asymptotic 
decomposition of the stress state starts with a solution describing the Saint-Venant edge effect in the theory of 
transversally isotropic inhomogeneous plates. In the case of transversal isotropy of the radially inhomogeneous material 
of the sphere, some boundary layer solutions decay very weakly, they can penetrate deep far from the conic sections and 
change the pattern of the stress-strain state. Asymptotic relations for displacements and stresses are derived, which 


provide calculating the three-dimensional stress-strain state of a radially inhomogeneous transversally isotropic sphere 
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of small thickness with any predetermined accuracy. It is shown that the root branching generates a countable set of 
new solutions for a transversely isotropic radially inhomogeneous sphere. 

Discussion and Conclusions. An asymptotic analysis of the stress-strain state of inhomogeneous shells, based 
on three-dimensional equations of elasticity theory, makes it possible to establish the limits of application of 
approximate theories. The identified behavior pattern of the solution far from the ends for different boundary conditions 
on the side surfaces can become the basis for creating refined applied theories for calculating the deformation of a 
radially inhomogeneous transversally isotropic spherical shell of small thickness. One of the applications of the 
asymptotic analysis performed can be the calculation of shells with thin coatings, in which, in this case, a radial 


inhomogeneity arises [13, 14]. 
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Introduction. Mechanical vibrations are widespread in the production processes. The drives of machines and 





mechanisms are mainly electromechanical, so mechanical reactive power is transformed into electrical reactive power 
of the network, impairing the quality of electricity. This explains the significance of considering the mechanical reactive 
power, and, as a consequence, the urgency of the presented study. The research objective is to detail the types of 
mechanical power under harmonic vibrations. 

Materials and Methods. The literature on the issues of dynamics, kinematics, vibrations, transformation of motion in 
oscillatory systems, etc., has been studied. Theoretical, mainly mathematical methods of research are used. 

Results. The powers developed under elastic deformations, forced harmonic vibrations of an inert body, and vibrations 
associated with gravitational influence, as well as reactive, active, full powers in the complex formulation, and 
mechanical powers in the vector representation are mathematically interpreted. 

Discussion and Conclusions. Under the mechanical harmonic vibrations, along with the sign-positive thermal power, 
sign-variable reactive powers develop, characterizing the reversibility of kinetic and potential energies. The total 
mechanical power satisfies the Pythagorean formula. The concept of mechanical reactive, active, and total powers 
generalizes the corresponding concepts of power from electrical engineering, and thus manifesting electromechanical 
dualism. 
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Introduction. Mechanical energy can be reversible (potential and kinetic), as well as irreversible (e.g., thermal 
in friction). The time derivative of the latter is taken as mechanical power. Note that due to the irreversibility of thermal 
energy, its derivative takes only positive values. At the same time, derivatives are obtained from both potential and 
kinetic energy. Of particular interest are harmonic vibrations [1—4], where derivatives (instantaneous powers) will be 
alternating functions, which fundamentally distinguishes them from thermal power. 

The analogue of kinetic energy in electrical engineering is the magnetic field energy of the inductor, the 
analogue of potential energy is the capacitor electric field energy, and the analogue of mechanical thermal energy is the 
thermal energy dissipated by the resistor. 

Mechanical vibrations are widespread in various production processes [5—8]. The drives of machines and 
mechanisms are mainly electromechanical [9—12], therefore, mechanical reactive power is transformed into electrical 
reactive power of the network, impairing the quality of electricity. In this regard, accounting for mechanical reactive 
power is of no small importance [14], and this determines the urgency of the presented study. 

Materials and Methods. Mechanical power under harmonic vibrations is considered. As a literary basis, 
domestic and foreign sources, which cover the issues on dynamics, kinematics, vibrations, transformation of motion in 


oscillatory systems, etc., have been studied. Theoretical (mainly mathematical) research methods are used. 
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Research Results 
The power developed by forced harmonic vibrations of an inert body. The body motion is described by the 
well-known expression: 
x=/sinat. 
Accordingly, the speed: 


v=x=lacosat =V, cosat. 


For harmonic quantity, the effective value is less than the amplitude in V2: 
Ve =—odl@ 


Formula for the force is: 
fi =m= —Ima@ sin at. (2) 
Formula for friction force is: 
Sf, = HX = LUa@cos at . (3) 
Resultant force is: 
f=f,+f, =—lma’ sin ot + ulacos ot = 
= lol we + n@ | ——“—— cos. ot - a — sin at J. 
(etm a (etna 
Denote: 
gy =arctg a, (4) 


Li 
With this in mind: 


T= len| ue +m or (cos pcos @t — sin gsin at) = lanl +m o° cos(at+¢@). 
FF. =loyl +m a. 


The effective value of the resultant force: 


ae ON a +m oO (5) 


V2 V2 


Obviously, 


Instantaneous resultant power: 


s= fv=lol iu +m’o cos(ot+@)lacos at = 
= 0,51" Vie +m ao [cos o+cos(2at + 9) | = 
= FV |cosg+cos(2at + v)|= 
= FV (cos p+ cos 2@t cos y — sin 2at sin g) = 
= FV cosp(1+cos 2@t)— FV sin psin 2at = p+ q,. (6) 
In electrical engineering, there is an expression similar to (6), with substitutions F — U V — I. The active 
power is determined from it: 
P=Ulcos@. 
Therefore, the active (thermal) mechanical power should also be defined as: 
P=FV cosg. (7) 
It is obvious that the harmonic force and velocity vibrate with a phase shift equal to g. 
From the above-mentioned formula of electrical engineering, the reactive power is determined: 
P=UlI sing. 
So, reactive (inertial) mechanical power should also be defined as: 
O,=FV sing. (8) 
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It follows from (6) that the active power is the half-period average of the instantaneous power, and the reactive 
power is the amplitude value. In electrical engineering, it 1s similar. 
Another generalization from electrical engineering is total mechanical power: 


S=FV =JO? +P”. (9) 


It is remarkable in that, on the one hand, it is described by the Pythagorean formula, and on the other hand, it is 
equal to the product of the effective values of harmonic quantities. 
With a view to (1), (5) and (8), 











loos we +mn’oa lo ma mo 
Q, = FV sn. g = — _—_ = = (10) 
2 /2 Ue +m" 2 
Here: 
f,v=—lma@ sin ot locos at =—0,51’ma” sin 2at =—F V sin 2@t =—O, sin 2at. (11) 


This corresponds to expressions (6) and (10). 
Keeping in view (1), (5) and (7), 








lanl we +mo lo Li ul o 
P= FV cosg = —___ = = = (12) 
V2 V2 fi +m'o® = 2 
In addition: 
f,,¥ = Hlocos at lacos at = 0, 5 ul @ (1+cos 2@t) = FV (1+ cos 2a@t) = P(1+cos 2a). (13) 


This corresponds to expressions (6) and (12). 
Having in view (9), (10) and (12), 
sa py love tm'o! lo _ Poti’ +m'ot 
V2 V2 2 
The power developed under elastic deformations. The expression for the force 1s: 
fF, =kx =klsin at. (14) 
Taking into account (3), the resultant force is: 


f=he +f, =hMsin ot + lacos wt = 


=I fk? + Po | 


__* gg, ot +S cos cot 
ke +o ke +u Oo 
Denote: 
k 
o = arctg —. 
Lo 


Hence, 


f=l\k’ +o (singsin at + cosgcosat)=l lk? + a" cos(at—¢@). 
It is obvious that: 
F=f + war. 


The effective value of the resultant force is equal to: 


F, Wk +o aS 
V2 V2 


Instantaneous resultant power is: 


s=fv=lik’ +u’@ cos(at—¢)lacos at = 
= 0,5? alk? + wa” | cos p+cos(2at-¢) |= 


= FV |cos g+cos(2at — )|= 


= FV (cos p+ cos 2a@t cosy + sin 2at sin g) = 
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= FV cos g(1+cos 2@t)+ FV sin gsin 2at = p+q,. (16) 


Bearing in mind (6), (7) and (12), the active mechanical power is equal to: 


Ifk? + -O 
P=FV cosg= ee ae = 


— uo 


o 
V2 V2 JP 4+p?o® = 2 


Taking into account (15), (1), (8) and (16), the mechanical reactive (elastic) power is equal to: 








QO, Spying Ae ee (17) 
2 2 (e+eo 2 
Here: 
fv =klsin ot locos at = 0,5kl’ asin 2at = FV sin2ot = QO, sin 2at. (18) 


This corresponds to expressions (16) and (17). 

Obviously, the total power is: 

i i a 2 222 
S=FV = lo? +P? ee ee a 

Power during vibrations associated with gravitational influence. When the suspended load is deflected by 

angle a , the moment occurs: 
M =megLa. 
Suppose 


a=a, sinat. 


A=a,ocosat =a,, ie cos at. 


q, =Ma=mgLa, sinat a, (a ot = 0,5ma;, /Lg* sin 2at. 


Then 


Instantaneous power has the form: 


Its amplitude and, accordingly, the reactive power of the gravitational effect is defined as: 
0, = 0,5ma; Le. 
Reactive, active and total power in a complex representation. In [15], it is shown that under an inert load, 
V. =V.e ae 
In this case, the instantaneous speed 1s equal to: 
v=V_cosat = ImV,. 

Formulas for the effective values are not fundamentally different: 

V =Vel"! F = Felt) | 
In electrical engineering, a feature of the complex representation is described in detail: when calculating the 


total power, one of the multiplied vectors must be conjugate. 
S= PV = Fel oI? = pei #0") — Yel = FV cos + jFVsing =P+jQ. 
This is an expression for an inert load. The elastic load differs in that the reactive power has the opposite sign: 
S= PY = Fel? Oye I? = FVel?-2-*) — Fe” = FV cos o— jFV sing = P+ jQ,. 
here: 
P=ReFV, Q=ImFV. 
Mechanical powers in vector representation. The complex representation is based on the concept of vectors 


rotating in the complex plane. The same principle can be implemented in a three-dimensional Cartesian basis. 
From (7)-(9), it follows: 


P=(F,V), O=|[F.V], S =(F.V) +[F.7]. 





Mechanics 


Zz} 


http://vestnik-donstu.ru 


i) 
ioe) 


Advanced Engineering Research 2022. V. 22, no. 1. P. 24-29. ISSN 2687-1653 





Mathematical abstraction with projections of rotating vectors has a specific material basis in the form of slot- 
and-crank mechanisms. 

Discussion and Conclusions. Mathematical methods were used to study the following powers: 

— under forced harmonic vibrations of an inert body, 

— under elastic deformations, 

— under vibrations associated with gravitational action, 

— reactive, active and total (in complex representation), 

— mechanical (in vector representation). 

It is shown that under mechanical harmonic vibrations, not only sign-positive thermal power develops, but also 
alternating reactive powers characterizing the reversibility of kinetic and potential energies. 

At the same time, the total mechanical power satisfies the Pythagorean formula. 

The idea of mechanical reactive, active and total powers is a generalization of the corresponding concepts of 


power from electrical engineering, and thus, electromechanical dualism manifests itself. 
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Introduction. The launch vehicle (LV) in flight and the dynamic components of loads from the impact of a trapezoidal 
wind gust are considered. It is proposed to determine the dynamic components of the force factors using analytical 
solutions for the structure points accelerations. The work objective is to create a technique for selecting the duration of 
the standard gust, under the influence of which maximum loads are provided in the sections of the LV structure. 
Materials and Methods. The launch vehicle is presented as an uneven beam. The description of its vibrations 1s reduced 
to a system of independent ordinary differential equations that determine the motion of an equivalent system of 
oscillators. The equation of oscillator vibrations under the action of a trapezoidal pulse load is solved by the overlay 
method, and it is reduced to the calculation of the Duhamel integral. It 1s proposed to get the parameters of an 
equivalent system of oscillators based on the results of the calculation of dynamic characteristics for a finite element LV 
model in the Nastran program. 

Results. Analytical relations for the LV structure point accelerations under the action of a trapezoidal wind gust are 
given. For the beam model, test calculations of accelerations were carried out according to the technique proposed in 
this paper. These data are compared to the results of finite element modeling. With the help of analytical solutions, 
dependences are constructed that determine the nature of the change in the magnitude of the bending moment for 
different sections of the launch vehicle when the duration of the wind gust varies. 

Discussion and Conclusions. The presented technique provides building an equivalent dynamic model of systems with 
a large number of degrees of freedom on the example of a LV and obtaining analytical solutions for accelerations of 
points of a mechanical system under trapezoidal external action. These solutions are applicable for the study of dynamic 
loads. The analysis results enable to select the duration of the wind gust, at which maximum loads are reached in the 
sections of the LV structure. Calculations based on the analytical solutions are very economical in terms of time spent. 
They can be used in design calculations for preliminary assessment of loading. 


Keywords: \aunch vehicle, beam model, oscillator, structure loads, Duhamel integral, dynamic loads, gust, analytical 
solutions, differential equation, vibrations, bending moment. 
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Introduction. During the operation of the launch vehicle (LV), loads occur in the elements of its design. We 
are talking about longitudinal and shear forces, bending and torsion moments. Data on these force factors are used for 
strength analysis in the design of new products, experimental development of the design [1], and adaptation of launch 
vehicles for a specific start-up [2]. Loads are divided into quasi-static and dynamic. The quasi-static ones arising in 
flight are calculated from the condition of dynamic equilibrium of the LV as a solid body, taking into account the 
permissible parameters of the ascent trajectory. 

Generalized beam models are usually used to calculate dynamic loads. Such loads are determined by the results 
of solving the equation of the elastic LV motion, which in general is a partial differential equation. Methods based on 
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the application of decomposition of the solution by the eigenvibration tones of the structure are able to provide a high 
speed of calculation in combination with sufficient accuracy of the results [3]. As shown in [4], using the method of 
decomposition by the forms of natural oscillations, it is possible to proceed to a system of independent ordinary 
differential equations. They describe: 

— centroidal motion of the LV, 

— rotation of the longitudinal axis of the LV relative to the center of mass, 

— transverse elastic modes of the LV in flight. 

The transition to independent equations describing elastic modes of the LV means that the distributed 
parameters of oscillators will be taken into account, each of which is a single-degree-of-freedom system. The motion of 
such an oscillator can be considered independently of the others, and a solution can be obtained for each of them using 
well-known methods of oscillation theory. 

To calculate the loads in flight, not the entire ascent trajectory is considered, but only some of its points, the so- 
called load cases, characterized by the extreme value of individual parameters affecting loading, or the maximum value 
of loads on individual structural elements. One of the most important cases of loading 1s the atmospheric flight of the 
LV [1, 5]. The influence of a turbulent atmosphere on the LV loading can be determined by statistical methods [6, 7], or 
within the framework of a conservative approach, when the maximum possible (with some level of probability) wind 
characteristics are taken into account. This paper discusses the second approach. A single specified wind gust is 
accepted as an external dynamic effect. The profile of the specified gust, which characterizes the change in wind speed 
over time, can be set in a trapezoidal [8], cosine, or sinusoidal form [9]. In this paper, we will consider the launch 
vehicle motion under the trapezoidal gust action. The duration of the specified gust is usually selected to be comparable 
with the period of the lowest transverse tone of the LV’ vibrations. At the same time, there are often demands for its 
variation to achieve maximum efforts in the LV sections [8, 10]. The complexity of calculations using standard finite 
element (FE) analysis programs is due to: 

— the need to vary the parameters of external action, 

— a large number of calculated cases, 

—a variety of design options and configurations of structures at the stage of design calculations [9]. 

This paper objective is to develop a methodology for selecting the duration of the specified gust using 
analytical solutions obtained for a simplified dynamic LV model presented as an equivalent system of oscillators. 

A semi-analytical approach using the Duhamel integral was successfully applied in [11] for hydroelastic 
analysis of ships. In [12] and a number of other works, the Duhamel integral is used as part of the problem solution of 
loading bridges with moving loads. In this paper, the Duhamel integral is used for analytical solutions to the LV 
reaction to the short-term impact of a wind gust in flight. 

Materials and Methods. At the stages of preliminary design, it is advisable to use flat design schemes for 
beam models. With the simplicity and speed of the solution, they enable to determine the motion variables and internal 
forces (with an accuracy acceptable for this stage of design) [13]. Imagine the LV in the form of an elastic beam with 
variable length mass and stiffness. Assume the usual assumptions for the resistance of materials, including the 
hypothesis for the smallness of elastic deformations. To determine the internal force factors in the LV section, we use 
the acceleration (overload) method [1, 4], which can be interpreted as a cross-section method adapted for dynamic 
calculation. In this case, internal forces are found from the conditions of static equilibrium of mentally cut off parts of 
the structure under the action of external distributed loads, supplemented by D'Alembert’s forces of inertia, and the 
desired internal forces. Quasi-static and dynamic values of force factors are determined separately based on pre- 
calculated accelerations, and then summed up [4]. 

This paper considers the issue of determining the dynamic loading of the LV in the transverse direction under 
the action of a wind gust, whose speed is directed perpendicular to the longitudinal axis of the LV. It is assumed that the 
loading in the longitudinal direction can be calculated independently. It is not considered in this paper. 

To determine dynamic accelerations, the LV is presented as a free elastic beam. Its motion is studied in the 
vicinity of the moment of time corresponding to the load case under consideration, and is described in deviations from 
the state of dynamic equilibrium in which the LV was before the wind gust, moving along the nominal (undisturbed) 
trajectory. Here, such parameters as the mass and moment of inertia of the LV, the motor power, and the angle of 
projection are assumed to be constant and equal to the characteristics of the considered point of the nominal trajectory. 
The perturbed motion of the elastic LV is investigated in a fixed coordinate system associated with the position that the 
LV occupied at the time of the calculation. The perturbed motion will be a combination of plane-parallel motion of the 
LV as a rigid body in the plane in which the dynamic load is applied, and elastic modes of the structure. We do not take 
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into account the control system response, i.e., for the stabilization machine, we assume a long delay time compared to 
the time of application of the dynamic load. In general, under the influence of a wind gust, together with elastic modes 
of the body, the LV starts to move as a rigid body. The transverse component of the aerodynamic force, which is 
considered proportional to it at a small angle of attack, changes its value due to the displacement of the structure in the 
direction of the wind gust and rotation relative to the forward flow vector. The projection of gravity on the transverse 
axis of the LV also changes. To properly account for these changes, the equations of LV motion must be integrated with 
the equations describing the logic of the automatic stabilization, which is impossible at the early stages of design. 
Taking into account the significant mass and the moment of inertia of the LV, we will consider small: 

— the angle of rotation of the LV as a rigid body for the time of calculation; 

— the rate of displacement of the LV center of mass in the direction of the wind. 

This enables to ignore the impact of the above changes. The angle-of-attack increment (and, consequently, the 
transverse component of the aerodynamic force) is considered to depend only on the magnitude of the wind gust speed 
given as a function of time. Thus, taking into account the accepted assumptions, the aerodynamic load in the transverse 
direction is a load distributed along the length of the beam with a time-dependent proportionality coefficient. The law of 
distribution of aerodynamic load along the length of the LV is determined experimentally and is considered to be 
known in advance. The law of change of the proportionality coefficient (angle of attack) from time to time is 
determined by the selection of the profile of the specified wind gust. 

The motion of the LV modeled as an elastic beam can be described using the well-known equation of forced 
transverse vibrations of the beam, written with the account 2 Voigt's hypothesis: 


m(x) 229 + (14+ h2)F [B@) 2] = a0), (1) 


where m(x) — mass per unit length; B(x) — bending stiffness; q(x,t) — distributed external load; h — friction 








factor. 
This equation should be supplemented with a boundary condition: the internal forces in the initial and final 
sections are zero. It means: 


2 (BG) 2) = 0, B(x) 2E9 = 0, mupyx=0,x =L, (2) 
where L — the LV length. 


In this paper, the aerodynamic force is taken as an external distributed load, which can be presented as a 
product of functions: 


q(x,t) = R(t)Y,(%), (3) 
where R(t) — the function that determines the temporal variability of the aerodynamic force and varies according to 
the trapezoidal law in accordance with the wind gust model adopted in this paper; Y,(x) — the function of the 


aerodynamic force distribution along the LV length. 

Let us consider free LV vibrations without taking into account friction forces (at q(x,t) = 0, h=0). We 
substitute the variable separation method y(x,t) = f(x):q(t). In this case, from equation (1) with boundary 
conditions (2), it is possible to arrive at an Rae differential equation 


*[B@) 2] - p?ma@f@) = 0, (4) 





with boundary conditions: 








S[B@) LO] = 0, Bo) LO =o, npux=0,x=1 (5) 


Solution (4) with conditions (5) is a classical Sturm-Liouville problem. Solving it, one can find a set of 
eigenforms f;(x) and eigenfrequencies p; of the beam in question Vj = 1,2,...). It is known’ that some solutions to 
system (4) correspond to zero natural frequencies. The forms corresponding to zero natural frequencies determine the 
translational motion of the LV as a rigid body together with the center of mass and rotation around the center of mass: 
f-1 = 1, fo = x — Xc, where x; — coordinate of the LV center of mass. 

It should be noted that the LV mass and stiffness characteristics most often have a piecewise constant pattern 
of distribution. In this case, the equations of form (1) and (4) should be written separately for each homogeneous beam 
section with boundary conditions at the junctions of the sections, as in the derivation of the ratios of the initial 
parameters method [1,14]. In the given paper, this entry is omitted, because the calculation of the dynamic 
characteristics (modal analysis) of structures is carried out numerically, using the finite element method. 

Let us imagine the forced vibrations of an elastic beam modeling the LV structure in the form of 
decomposition according to the natural modes. Assume that the beam stiffness axis passes through its center of mass. 
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To move the points of the LV axis, we write: 

y(x,t) = ye(t) + 9) — xc) + D8 HR) “aj ©, (6) 
where y(t) — displacement of the beam center of mass; 9(t) — angle of rotation of the axis of the undeformed beam; 
f; — beam eigenmode corresponding to tone numbered j; q,;(t) — generalized coordinate corresponding to tone 
numbered j; N — number of elastic tones taken into account. 

After substituting (6) into (1) and applying the Bubnov-Galerkin procedure, we can arrive at an ordinary 
differential system with constant coefficients: 
MYc(t) = Qy, 
19(t) = Qo. (7) 
m,(q;(t) + 2njq,(t) + p74q;(t)) = Q(t) G = 1,2,...,N) 
Here, m — LV mass; J — moment of inertia relative to the axis passing through the LV center of mass perpendicular to 
the plane of rotation; m; — reduced (generalized) mass for the j-th vibration tone and determined from the formula: 


L 
mj = J, m(x)f/ (x) dx. 
The generalized forces in expression (7), taking into account (3), are defined as follows: 


i 
Q, = R(t) i Y, (x) dx = RU); (8) 
L 

Qs = R(t) [. Yq (x) (« — xX¢)dx = R(t)Mg, (9) 

i 
Q; = R(t) fe Y,(x) fjdx = R(t)Qo;. (10) 
Here, Q, — maximum value of the main transverse aerodynamic load vector; M, — maximum value of the main 
transverse aerodynamic moment reduced to the LV center of mass; Q); —- maximum value of the generalized force 


corresponding to the generalized coordinate qj. 

The first two equations in (7) define the law of change of accelerations of the LV points in the process of 
translational and rotational motion of the LV as a rigid body. The last equation in (7) defines the law of motion of an 
equivalent system of oscillators. 

Consider the motion of one oscillator under the action of a trapezoidal external load, which: 

— increases from zero to Qp over time 6, 

— maintains a constant value over time 0, 

— drops to zero over time 6. 

For convenience, we omit the indices characterizing the tone number. Then we will rewrite the differential 
equation of the oscillator motion taking into account (10) in the form: 


G+ 2ng + p?q = =“ RO), (11) 
where m — oscillator mass; p — angular frequency of natural oscillations, expressed in radians per second; n — 
damping coefficient (determines the oscillator damping); Q(t) = Q)R(t) — law of variation of external load. 

We represent function R(t) as a set of four linear functions: 

Re) =o. 2 R(t) = —8"'(t — 9), 
R3(t) = —6°-*(t — 6-8), R,(t) = 5 *(t — 26 — 9). 

Accordingly, external load Q(t) is a combination of four linear loads Q;(t) = Q)R;(t), (i = 1, 2,3,4). Load 
Q,(t) is applied from moment t = 0; Q,(t) — from moment t = 6; Q3(t) — from moment t = 6 + 6; Q,(t) — from 
moment t = 8 + 26. 

To determine the system response to external actions, we divide the entire duration of the load into four 
intervals (Fig. 1). 


Q(t),N 





' Ist section | 2nd section : 3rd section; 4th section 


Fig. 1. Diagram of the external trapezoidal load 
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In accordance with the superimposition method’, let us represent the response of the linear system under 
consideration to an external action in the form of the Duhamel’s integral as the sum of responses to a set of 


independently applied elementary impulses: 
t 


a(t) = | P(e- 8) (0)ab, 
0 
Here, P(t — 8) — the shifted time law of variation of external action, and Y(9) characterizes the system 


response to a unit impulse input. The response of mechanical system q,(t) to a linearly increasing load P(t) = kt can 


be expressed through the system response to unit impulse Y(9) and to a suddenly applied single load Y, (9): 
t 


q(t) =k | Y,(9)d8, Y,(t) = | Y(9)d9. 
0 


0 


The response of a one-degree-of-freedom mechanical system with damping coefficient n to a unit pulse will 


have the form’: 


1 
Y(9) = —————-e™”* sin ( p2 — n29), 


We introduce a notation for the frequency of damped vibrations p, = ,/p* — n2. Let us calculate the system 


response to a suddenly applied unit load: 
t 








1 n 
Viiv) = | Y(9)d9 = E —e (cos( t) + —sin( ‘))} 
1 mp? P41 D, P41 
The response to a linearly increasing load: 
_ t _  k (2n =i at (an? —nt o: 
q(t) =k in Y,(9)d9 = aa (s [—1 +e cos(p,t)] Les ( 2 1) e~™ sin(p,t) + th (12) 


We accept the notations: A,;= t — 6, A,= A, — 98, A3= A, — 4. 


We introduce a function containing the harmonic terms of solution (12): 
n he re yh 
d(9) = 2= (cos(p, -0) +—sin(p, -9)) e ™ ——sin(p,:8)e-™. 
p P41 P4 


We denote by q,, the motion under the action of a force statically applied to the system Q) = g,,mp” and take 
into account that k = Q,)/6. The total movement of the oscillator under the action of a combination of loads Q;(t) at 
every time point will be the sum of the corresponding solutions (12). The system response described by equation (11) to 


the external impact of the trapezoidal profile will have the form: 


daS "ld = 2n/p> +t] at0<t <6, 
oe derS *[d(t) — d(A,) + §] at8<t<6+6, 13) 
ded 1 [d(t) — d(A,) — d(A,) + 2n/p?—A,] atO9+5<t< 0+ 26, 
der *(d(t) — d(A,) — d(A2) + d(A3)) att > 0 + 26. 


The law of variation of the oscillator accelerations can be obtained from double differentiation in time of 


d2 


aa (d (9)). We differentiate and record the result taking into account 


expression (13). Let us introduce function g(9) = 


the number of the oscillation tone: 


Nn; Nn; 1 
g;(8) = G = pt ;) (24 (costo, -9)+ —-sin(p, ; )) = —sin(p, ; )) — 
Dj P41 j P41 j 
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=on; (22 (py, sin(p, ; -9) ‘Pad; cos(p, ; -9)) = cos(p, ; -8))| ew. (14) 


We rewrite the expression for the frequency of damped natural vibrations: 


Pij= Ip; ie (15) 


To determine the accelerations of the LV points, we differentiate twice (6): 


W(t, x) = ¥e(t) + OO) — xe) + VL HOG (O. (16) 
Let us express the accelerations of the generalized coordinates from the first two equations of system (7) taking 


into account (8) and (9). We represent the accelerations of the points of the LV axis that it acquires through moving as a 


rigid body: 
a(t) = ¥c(t) + 8(E) (x — xc) = R(t)(Qa/m+t M,/1). (17) 
In addition, we take into account the law of variation of function R(t) from time: 
5 1t atO <t <6, 
R(t) = 1 até6<t<68+46, 
—§ 4(t — 26-6) at0+6<t<0+4+26, 
0 att >0+ 26. 
We get: 
5-'t(Q,/m+ M,/1) ac 0=t.< 5, 
ioe Q,/m+M,/1 atOé<t<0+46, (18) 
—5-1(t — 26-6)(Q,/m+M,/I) w8+6<5t<0+ 26, 
0 att. 0+ 20. 


For accelerations of generalized coordinates corresponding to the elastic vibration tones, as a result of double 
differentiation of expression (13), taking into account (14), we obtain: 
5 der 59; (t) at0<t <6, 
8 qer j (gj(t) — 9;(Ay)) até6<t<6+6, 
8 der j(gj(t) — gj(Ar) — gj(A2)) at +5 <b <0 +26, 
8 der j (gj(t) — 9j(Ar) — gj(A2) + gj(A3)) att > 0+ 28 


Taking into account (16) and (17), the law of variation in the accelerations of the elastic LV axis points under 


qj(t) = (19) 


the action of a trapezoidal wind gust will have the form: 
VE, x) = a(t) + Yj ffO)GO, (20) 
where a(t) is determined from expression (18), g;(t) — from expression (19). 
While investigating the LV elastic modes, damping is traditionally taken on the basis of the data obtained from 
the results of full-scale dynamic tests and presented in the form of values of logarithmic decrements D;. Then, the 
coefficient determining the damping parameter and included in expressions (14) and (15) can be calculated from 


formula: 
The change in generalized coordinate q; included in formula (19) under the action of a statically applied 
generalized force Qo; is determined from expression: 


4 Qoj 4 fy Yai) ax 
Pp Mj = PF Jo mx) F? (x) dx 


cr J ~ (22) 


Taking into account (14), (15), (18), (19), (21), (22), formula (20) is an analytical expression that defines the 
functions of the acceleration change over time for the points of the LV axis under the action of an external transverse 
aerodynamic force varying according to the trapezoidal law, taking into account the impact of dissipative forces. 


Knowing the law of acceleration variation, it 1s possible to determine the dynamic, and then the total structure 
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loads acting in the LV sections using the known methods [4]. We calculate the bending moment due to inertial forces 
from the elastic modes of the LV structure: 
M(x,t) = — Vj Mj (1) GO. (23) 
Here, M;(x) — function of distributing a unit (with acceleration q;(t), equal to one) bending moment for the j-th 
vibration tone along the LV length. It can be found from formula: 
M,(x) = fr fo m(x)f; (x)dx dx. (24) 

As noted above, for the modal analysis in this work, the FE method was used. This approach is due to the fact 
that in practice, the dynamic LV model has rather complicated structure. It includes substructures, and their own 
dynamics cannot be neglected. Substructures can be attached to the LV body in one section, or be located parallel to the 
longitudinal axis of the LV and have several attachment points. In this case, the calculation of dynamic characteristics in 
a continuum setting is a complex mathematical problem. In addition, dynamic models of individual substructures are 
presented by development companies in a condensed (matrix) form in the Nastran format. For this reason, it will be 
optimal to use the Nastran engineering analysis software package to calculate the dynamic characteristics of the 
structure. 

However, the use of standard FE analysis programs for the calculation of dynamic loading is fraught with 
certain difficulties. These include the need to pre-construct an equivalent model of the external aerodynamic load 
suitable for use in the FE analysis program [15], which in itself is quite difficult. In addition, there are difficulties 
associated with processing the calculation results. The use of the postprocessor functionality for analyzing the results is 
extremely time-consuming and requires a large number of manual operations. Another way involves the application of 
additional software for processing large array of numeric data”. In this paper, an approach is proposed in which the FE 
analysis program is used only for modal analysis. In this case, dynamic loading is calculated using specially developed 
software that allows you to vary external loads and automatically process the calculation results. 

For a complex LV design, equations (7) retain their form [3]. The standard output information of the Nastran 
program can be the basis for obtaining the parameters of an equivalent system of oscillators, LV mass and moment of 
inertia, as well as for calculating the generalized forces included in the third equation of system (7). To form the left 
side of the third equation in (7), the values of natural frequencies p, (Radians) and generalized masses m, (Generalized 
mass) are required. To determine the generalized forces in the right part of the third equation in (7), the eigenmode 
functions f;(x) (Eigenvector) are needed. When calculating dynamic inertial loads according to (23), instead of (24), it 
is more convenient to use unit inertial loads (forces and moments), which are output by the Nastran program after the 
standard application of forces when calculating eigenforms and frequencies. Unit inertial loads are output separately for 
each vibration tone and multiplied by the square of the natural frequency, which should be taken into account for their 
correct use. 

Research Results. To carry out test calculations, a medium class tandem launcher is considered. Figure 2 
shows the type of functions of linear bending stiffness B(x) and mass m(x), as well as the distribution of concentrated 


MASSES Megcp(X) = Meocp A(X — x,) along the length of the considered LV (the Dirac delta function is denoted by A). 
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m(x), kg/m 


Mecocp (x) ) kg 


kg-N 
B(x), =— 


c) 


Fig. 2. Distribution of mass and stiffness characteristics along the LV length: 
linear mass (a), concentrated mass (b), bending stiffness (c) 


To check the obtained analytical solutions under the impact of a trapezoidal external load, the dynamic LV 
accelerations are calculated according to formula (20), taking into account (14), (15), (18), (19), (21), (22). At the same 
time, 5 elastic tones of natural transverse LV vibrations are considered. Characteristics p;, m, and f;(x) were obtained 
from the calculation results in the MSC Nastran software package using the solution sequence for modal analysis of 
natural vibrations (SOL 103). The LV dynamic finite element model is presented as a set of beam elements with various 
inertial and stiffness characteristics. The following elements are elastically or rigidly attached to the beam elements: 

— elements describing the inertial properties of devices, aggregates, parts of the block structure; 

— condensed models of individual blocks presented in the digital matrix form. 

The finite element LV model includes about 1000 components. The elements simulating restraints were not 
used to preserve the ability of the LV to move as a rigid body. 

According to the same finite element LV model in the MSC Nastran software package, the accelerations are 
calculated through the SOL 119 solution sequence used for modal transient analysis. At the same time, all tones of 
natural vibrations in the range up to 100 Hz were taken into account in the modal decomposition. The aerodynamic load 
is represented by transverse linear loads distributed over all beam elements simulating the LV structure. In addition, the 
trapezoidal law of the aerodynamic load variation over time 1s given. 

Figure 3 shows the results of a comparative analysis of accelerations obtained using two different approaches 
A(t) = V(t,xX9) of a certain point of the LV axis with coordinate x = x9. It can be seen that the following two 


solutions agree well: 
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— the solution obtained for accelerations using a simplified LV model based on the analytical relations given 


in this paper; 
— numerical solution obtained from the complete finite element LV model. 
a,(t), 30.0 
2 
m/s" 20.0 


10.0 
0.0 + 
-10.0 
-20.0 
-30.0 + 





Fig. 3. Comparison of transverse accelerations obtained through two different 
— finite element modeling, ® °° °° — analytical solution for a simplified model 





approaches: 

The peak value of the bending moment selected from the time process determines the level of equivalent forces 
taken to carry out the strength calculation, and acts as a variable parameter when changing the parameters of external 
action in the transverse direction [16]. 

Dependences are obtained based on the values of dynamic accelerations (Fig. 4). They show how the change in 
peak bending moment M or different LV sections are related to parameter 8, which characterizes the gust duration. In 
Figure 4, the bending moment values are presented in the form of dimensionless quantities M*. They are calculated 
through dividing the dimensional bending moment by the maximum value for a given section (e.g., for section x = 0.3L 


maximum value M,,,,, = 4.6: 10° N-m), found when value 0 varies over the entire range under consideration. 


M" 1.0 





a) 








c) 


Malykhina O. I. Analytical solution to approximate equations of the launch vehicle motion under the gust action 











Fig. 4. Dependence of the value of dimensionless bending moment M* on the duration of wind gust 0*, expressed in fractions of 
period T, for various LV sections: zone 1 (x = 0.0 — 0.2L) (a); zone 2 (x = 0.2L — 0.4L) (6); 
zone 3 (x = 0.4L — 0.55L) (c); zone 4 (x = 0.55L — 0.75L) (d); 
zone 5 (x = 0.75L — 0.9L) (e); zone 6 (x = 0.9L — L) (f) 


So, on the graphs, the duration of the gust action is represented by dimensionless value 0° obtained through 
dividing parameter @ by period T, of the first tone of the LV vibrations. Each line of the graph corresponds to one 
section of the LV. All LV sections are grouped according to the nature of function M*(6*) and are shown on various 


graphs, and the LV length is appropriately divided into zones (Fig. 5). 






— 


Fig. 5. LV cross-section zones 





It can be seen from Figure 4 that for the first zone, the maximum bending moment is reached already at 0 
values at least 15% of period T, of the first vibration tone. For the second and third zones — slightly more than half of 
period T,. A further increase in 0 values does not affect the magnitudes of the maximum values of the bending moment. 
For zones 4-6, the maximum value of the bending moment turns out to be local and is within the zone of 8 values close 
to the value of half of period T, of the first tone of the LV vibrations. The results obtained fully correspond to the results 


of finite element modeling carried out earlier®. 


° Malykhina OI, Glugovskii MS. Analiz vliyaniya profilya poryva vetra na velichinu korpusnykh nagruzok rakety-nositelya v poletnykh sluchayakh 
nagruzhentya. In: Proc. VI All-Russian Sci.-Tech. Conf. Samara. 2019;1:133-138. (In Russ.) 
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Thus, to obtain the maximum values of the bending moment in the sections of the considered LV, a dynamic 
analysis of the behavior of the structure with external action in the form of a wind gust of is required. Its duration is 
determined by parameter 9, close in value to half of the period of the first tone of the LV vibrations. 

Discussion and Conclusions. Using the superposition method, analytical solutions are obtained that describe 
the motion of a one-degree-of-freedom system, which is affected by the friction force and an external force varying 
according to the trapezoidal law. The method of application of the obtained analytical solutions for systems with many 
degrees of freedom is given. A good coincidence of two types of solutions is shown: 

— analytical one, for accelerations of LV points found from a simplified model; 

— numerical data obtained from the complete finite element LV model. 

It is shown that analytical solutions can be used to analyze dynamic force factors to select the duration of a 
wind gust, under the impact of which maximum loads are achieved in the sections of the LV structure. Similarly, it is 
possible to analyze the overloads that are achieved in the LV sections (e.g., at the installation points of measurement 
systems). 

In addition, the proposed methodology provides building a full cycle of the load analysis pre-calculation in the 
case when an analytical representation of the external dynamic load is possible. The load analysis based on analytical 
solutions 1s very economical in terms of calculation time, and it can be a remarkable alternative to finite element 
modeling at the design stage, when a large number of combinations of external loads and configurations of the design 
under development are studied. Finite element analysis of the detailed model in this case can be used as a refine final 


calculation. 
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Introduction. To search for effective methods of finishing and stripping of long-length parts, it is advisable to create a 
rotary-screw process system of a pass-through type. Its main working element is a screw rotor, which is a combination 
of flat elements of various shapes and sizes, multidirectional with respect to the helical lines around the perimeter. The 
purpose of this study is to justify the intensity of the machining process in devices equipped with a screw rotor. 
Materials and Methods. Metal removal is accepted as the main parameter determining the intensity of the machining 
process in screw rotors. Within the framework of the presented study, processing was performed on an experimental 
rotary-screw installation. The processing medium consisted of molded abrasive pellets of the PT 10x10 brand. The 
research conditions were as follows: loading volume (without part) — 60 %; rotor speed — 50 rpm; processing time — 
30, 60, 90 min; rotor axis tilt angles — 0° and 5°. The influence of treatment modes and conditions on the process 
intensity was considered on plate samples with dimensions of 80x10x1 mm made of aluminum alloy DI6T. To 
determine metal removal, the samples were weighed on Ohaus AX223 analytical balance before and after processing. 
Results. The patterns of metal removal from samples at different processing times, the location of samples in the 
working area of a screw rotor with different angles of inclination, and its speed are presented. 

Discussion and Conclusions. The regularities established in the course of research indicate the efficiency of rotary- 
screw process systems for solving problems on the finishing treatment of long-length parts. The main factors that 
control such processing in devices with a screw rotor and affect its intensity (and, as a consequence, the process 


performance) are the screw rotor speed, the shape of the perimeter, and the axis angle. 
Keywords: finishing and cleaning, screw rotor, long parts, processing medium, process intensity, metal removal. 
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Introduction. Finding efficient ways of finishing and cleaning long parts [1-6] predetermined the creation of a 





rotary-screw through-flow process system to solve this problem (Fig. 1) [7-15] "”. 


' Babichev AP, Babichev IA, Serga GV. Device for vibration processing of long parts. RF Patent no. 2228252, 2004. (In Russ.) 
* Lebedev VA, Al-Obaidi LMR, Koval NS, et al. Unit for finishing and cleaning of long parts. RF Patent no. 2750922, 2021. (In Russ.) 
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Fig. 1. Process flowchart of finishing and cleaning long parts in a screw rotor: 
1 — mounting frame, 2 —lugs, 3 — adjusting screws, 4 — engine, 5 — sump tank, 6 — drive shafts, 7 — roller bearings, 
8 — shells, 9 — rollers, 10 — rims, 11 — screw rotor, 12 — gasket, 13 — slotted gasket, 14 — piping system, 15 — part, 
16 — processing medium, 17 — container 
The main working element of this type of process system is a screw rotor. Its perimeter is made in the form of 
a combination of flat elements of various shapes and sizes, multidirectional with respect to the helical lines on their 
outer surface. The processing granular medium is located inside the rotating screw rotor. Being carried away by flat 
elements, at a certain moment, under the influence of gravity, it avalanches downwards. At the same time, the upper 
layers of the loading masses rotate around their own axes, drawing in nearby layers of particles of the loading masses. 
As a result, the particles of the loading masses roll relative to each other. In addition, when an avalanche rolls down the 
uneven surface of the underlying layers, small force strikes and sliding of the particles of the loading masses occur. 
Thus, the part is processed by abrasion, scratching, and, to a small extent, by blows. The technological effect of 
processing depends on the dimensions of the screw rotor, the particle masses of the processing medium, the filling 
factor of the working chamber, the speed of rotation, the geometry of the screw rotor, and a number of other factors. 
The objective of the presented study is to justify the intensity of the machining process in devices equipped with a 
screw rotor. 
Materials and Methods. Metal removal is accepted as the main parameter determining the intensity of the 
machining process in screw rotors. 
For processing, an experimental rotary-screw installation, made on the basis of a 1K625 model screw-cutting 


machine, was used (Fig. 2). 





Fig. 2. General view of the experimental rotary-screw installation 
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The processing medium consisted of molded abrasive pellets of the PT 10x10 brand. The research conditions were as 
follows: 

— loading volume (without parts) — 60 %, 

— rotor speed — 50 rpm, 

— processing time — 30, 60, 90 min, 

— rotor axis tilt angles — 0° and 5°. 

The effect of processing modes and conditions on the process intensity was considered on plate samples with 
dimensions of 80101 mm made of aluminum alloy D16T. Blanks in the form of long-length rolled products are made 
of this material, requiring finishing and cleaning before going into production. 

In the working area of the experimental screw rotor, three sections 4, B, C (Fig. 3), remote from the left end of 


the screw rotor at distances I) l, lz, were determined along the length of the working area. 




















Fig. 3. Diagram of the installation of samples on the mandrel in the rotor: 
1 — samples; 2 — mandrel; 3 — lathe toolholders 


The samples were fixed on a long four-sided hollow rod and installed according to the scheme shown in 
Figure 3. Distances l) 1, 1, correspond to the position of the sections of the rotor working area. To assess the nature of 
the impact of the processing medium on the various surfaces of the mandrel, 4 samples were processed in each position 
(section). They were fixed on the edges of the mandrel in such a way that samples | and 2 were in the zone of active 


exposure to the processing medium (with a sliding layer), and samples 3 and 4 were in the passive zone (Fig. 4). 








Fig. 4. Layout of samples on the mandrel in cross section: 
1, 2, 3, 4 — samples; 5 — screw rotor; 6 — processing medium 


The mandrel with samples was installed in the lathe toolholders and inserted into the working area of the screw 


rotor for processing. 
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To determine the mass removal of metal, samples were weighed on Ohaus AX223 analytical balance before 
and after processing. This model is designed for static measurements of the mass of substances and materials accurate to 
0.001 g. Metal removal was controlled after thorough washing and drying of the samples. 

The rate of metal removal from a unit surface area equal to the area of the square packing granules of the 


processing medium was determined from the ratio: 


Vya = St 





, 2/S, 


where AQ — metal removal from the surface, g; t — processing time, s; S — area of the processed surface of the 
samples, mm’; Syn — square packing granules of the processing medium, mm’. 

The processing duration was recorded using a time relay (timer) of VL-45UHL41 model. 

The experimental results were determined as the arithmetic mean difference of the corresponding measurement 
values of the samples or batch of samples before and after processing. 

Research Results. Figures 5—7 show the results of experimental studies on the intensity of the processing in a 


screw rotor. 
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Fig. 5. Change in material removal from D16T samples depending on the 
processing time and the location of the samples on the mandrel: 
a) with an angle of inclination of 0°; b) with an angle of inclination of 5° 
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Fig. 6. Change in material removal from D16T samples depending on 
their location in the working area of the screw rotor (in 90 min): 
a) with an angle of inclination of 0°; 5) with an angle of inclination of 5° 
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Fig. 7. Change in specific metal removal from D16T samples depending 
on their location on the mandrel (in 90 min): 
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Puc. 8. Change in metal removal from D16T samples depending on the loading of the screw rotor working area by processing 
medium and the rotation speed (in 30 min) 


Discussion and Conclusions. The research results allow us to make a number of statements. 

1. The form of the curve describing the change in the amount of metal removal from the surface on the 
processing time is close to linear (Fig. 5). This does not depend on the location of the samples in the working area of the 
screw rotor, its design, and the angle of inclination relative to the axis of rotation. 

2. The kinetics of particle motion and, as a consequence, the intensity of metal removal in various sections of 
the working area of the screw rotor are determined by the angle of inclination of the plates that create the shape of the 
working sections of the rotor. With their identical shape (Fig. 3), almost the same processing intensity is provided. The 
difference in the amount of removal in three sections was 3 %. 

3. Processing intensity (Fig. 5—7) depends on the position of the surface of the samples in relation to the 
direction of motion of the granules of the sliding layer (Fig. 4). Two cases of metal removal were compared; first— 
from samples | and 2 installed in the active zone of the processing medium; second — from samples 3 and 4 installed 
in a passive zone characterized by weak energy-force interaction with the treated surface. In the first case, the metal 
consumption is twice as high. 

4. Changing the angle of inclination of the screw rotor axis relative to the axis of its rotation increases the 
intensity of sample processing in the middle section (B—B) by 1.2 times and reduces it in the working areas of the end 
sections (Fig. 5—7). This is due to a change in the kinematics of the medium motion in the working area of the screw 
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rotor. With every full turn, the processing medium tends to move to the middle part of the rotor working area 
(section B—B). This provides: 

— constancy of the media volume pressure on the samples; 

— activation of shock-pulse action of the processing medium granules on the surface of samples. 

In addition, when the angle of inclination of the screw rotor axis changes, the axial pressure of the processing 
medium on the end walls decreases. This is very important for pass-through type installations, because this prevents the 
granules from spilling out of the holes designed for the passage of long parts through the working area. 

5. The rotation speed of the screw rotor is an important technology factor that determines the processing 
intensity in a rotary-screw installation. It can be seen from Figure 8 that with different loading degrees, the amount of 
material removal from all samples increases with an increase in the rotation speed of the screw rotor up to 50 rpm. 
When the speed increases to 70 rpm, the indicator decreases. This is due to the increased centrifugal forces. Their action 
causes a certain volume of particles of the processing medium to remain permanently at the rotor walls. As a result, the 
number of sliding layers of the medium is reduced, and, consequently, the intensity of the impact on the surface of the 
samples 1s reduced. 

The regularities and technological capabilities of the processing in a screw rotor established in the course of 
research indicate the acceptability and efficiency of the rotary-screw process systems for solving problems related to the 
finishing processing of long-length parts. 

In the conducted study, the key factors considered were the rotation speed of the screw rotor, the angle of 
inclination of its axis, and the shape of the perimeter. They provide finishing control in devices with a screw rotor, 
affect its intensity, and hence, the process performance. 


References 

1. Babichev AP, Motrenko PD, Gillespi LK, et al. Primenenie vibratsionnykh tekhnologii na operatstyakh 
otdelochno-zachistnoi obrabotki detalei. Rostov-on-Don: DSTU Publ. Centre; 2010. 285 p. (In Russ.) 

2. Babichev AP, Babichev IA. Osnovy vibratsionnoi tekhnologii. 2nd revised and enlarged ed. Rostov-on-Don: 
DSTU Publ. Centre; 2008. 693 p. (In Russ.) 

3. Butenko VI. Finishnaya obrabotka poverkhnostei detalei: sposoby, ustroistva, instrumenty. Rostov-on-Don: 
DSTU Publ. Centre; 2016. 219 p. (In Russ.) 

4. Tamarkin MA, Tishchenko EE, Shvedova AS. Optimization of Dynamic Surface Plastic Deformation in 
Machining. Russian Engineering Research. 2018;38:726—727. http://dx.do1.org/10.3103/S1068798X 18090277 

5. Tamarkin MA, Tishchenko EE. Osnovy optimizatsii protsessov obrabotki detale1 svobodnym abrazivom. 
Saarboniken: Lambert Academic Publishing; 2015. 140 p. (In Russ.) 

6. Li Xin, Peng Gaoliang, Li Zhe. Prediction of seal wear with thermal structural coupled finite element 
method. Finite Elements in Analysis and Design. 2014;83:10—21. http://dx.doi.org/10.1016/j.finel.2014.01.001 

7. Taratuta VD, Belokur KA, Serga GV. Rotary screw technology systems for combine harvesters. Proceedings 
of the Kuban State Agrarian University. 2015;57:197—206. 

8. Pesin MV. Improving the Reliability of Threaded Pipe Joints. Russian Engineering Research. 2012;32:210— 


212. http://dx.do1.org/10.3103/S 1068798X 12020232 
9. Macdonald KA, Bjune JV. Failure analysis of drillstrings. Engineering Failure Analysis. 2007;14:1641-— 


1666. http://dx.do1.org/10.1016/j.engfailanal.2006. 11.073 

10. Xiao-Hua Zhu, Yu Wang, Hua Tong. The parameter sensibility analysis for fishing box tap based on the 
overall process of elastoplasticity in oil and gas wells. Mathematical and Computer Modeling. 2013;58:1540—1547. 
http://dx.doi.org/10.1016/j.mcm.2013.06.004 

11. Fares Y, Chaussumier M, Daidié A, et al. Determining the life cycle of bolts using a local approach and the 
Dang Van _ criterion. Fatigue & Fractureof Engineering Materials & Structures. 2006;29:588—596. 
http://dx.doi.org/10.1111/j.1460—2695.2006.01029.x 

12. Serga GV, Seryy DG, Marchenko AYu. Investigation of physical phenomena occurred in contact area of 
bulk particles at their motion in screw drums by methods of similarity theory, engineering and computer graphics. 
Vestnik Bryanskogo gosudarstvennogo tekhnicheskogo universiteta. 2019;79:20-28. 
http://dx.do1.org/10.30987/article 5d10851f18f085.56011612 

13. Serga GV, Lebedev VA, Belokur KA, et al. Increased productivity of technological systems for finishing 
and strengthening parts based on using screw rotors. Strengthening Technologies and Coatings. 2016;136:16-19. 

14. Lebedev VA, Serga GV, Chaava MM, et al. The study of fine-cleaning treatment for removal of burrs in 
screw rotor. IOP Conference Series: Materials Science and Engineering this link is disabled. 2021;1029:012001. 


http://dx.do1.org/10.1088/1757-899X/1029/1/012001 


Lebedey V. A., et al. The finishing and cleaning of long parts in screw rotors 





15. Lebedev VA, Serga GV, Chaava MM. Method of Calculating the Machines Drive with Screw Working 
Bodies Mounted from Tetrahedral Hollows. ICIE 2021: Proceedings of the 6th International Conference on Industrial 
Engineering (ICIE 2020). 2021. P. 557-563. http://dx.do1.org/10.1007/978—3—030—54814—8 64 


Received 29.12.2021 
Revised 24.01.2022 
Accepted 25.01.2022 


About the Authors: 
Lebedev, Valery A., professor of the Mechanical Engineering Department, Don State Technical University 
(1, Gagarin sq., Rostov-on-Don, 344003, RF), Cand.Sci. (Eng.), professor, ResearcherID, ScopusID, ORCID, 


va.lebidev@yandex.ru 


Al-Obaidi Luay Mohammed Rajab, postgraduate student of the Mechanical Engineering Department, Don 
State Technical University (1, Gagarin sq., Rostov-on-Don, 344003, RF), lecturer, Middle Technical University, 
Institute of Technology-Baghdad (Muasker Al Rashid Street, Baghdad, Iraq) Scopus, ORCID, luayrajab@yahoo.com, 


luay_mohammed@mtu.edu.ig 


Koval, Nikolay S., associate professor of the Tool Engineering and Biomedical Engineering Department, Don 
State Technical University (1, Gagarin sq., Rostov-on-Don, 344003, RF), Cand.Sci. (Eng.), associate professor, 


Researcher ID, ScopusID, ORCID, koval-nc@mail.ru 


Claimed contributorship 

V.A. Lebedev: academic advising; basic concept formulation. N.S. Koval: research objectives and tasks 
setting; computational analysis; text preparation; formulation of conclusions. L.M.R. Al-Obaidi: analysis of the research 
results; the text revision; correction of the conclusions. 


All authors have read and approved the final manuscript. 


Machine building and machine science 


Bs 
\O 


http://vestnik-donstu.ru 


Advanced Engineering Research 2022. V. 22, no. 1. P. 50-56. ISSN 2687-1653 





MACHINE BUILDING AND MACHINE SCIENCE 





® Check for updates 


UDC 621.9 + 621.2.082.18 Original article 
https://doi.org/10.23947/2687-1653-2022-22-1-50-56 








Tribotechnical properties of experimental hard alloys with modified cobalt 
binder 


re 
[al 
Evgeniy V. Fominov¥> ><], Constantine G. Shuchev @® , Muharbiy M. Aliev © 


Don State Technical University (Rostov-on-Don, Russian Federation) 


bX] fominoff83 @mail.ru 


Introduction. This paper discusses tribomechanical characteristics of experimental hard alloys with a modified cobalt 
binder under friction without lubrication on hard-to-cut materials — stainless steel and titanium alloy. The research 
objective is to evaluate the process of friction interaction for each friction pair according to a number of parameters, and 
to determine the optimal combinations of “experimental hard alloy — structural material” on the basis of the established 
tribological indicators. 

Materials and Methods. Tribological tests of hard alloys were carried out using a cylinder-to-disc friction scheme for 
different sliding speeds and temperatures under constant load without the use of lubricants. Comparison of the friction 
interaction process was carried out by the frictional force, volumetric wear and roughness of the friction tracks on the 
counterbody. Stainless steel 12H18N9T and titanium alloy BT3-1 were used as counterbody materials. The resistance of 
experimental compositions to the abrasive type of wear was determined through measuring the surface dynamic 
microhardness on a scanning nanohardness tester by analyzing the thickness of the scratches caused by the indenter. 
Results. According to the results of surface microindentation, the experimental alloys 2.22 (binder 5.65% Co + 1.8% Mo 
+ 0.6% Ti) and 2.23 (binder 5.1% Co + 2.7% Mo + 0.61 % Ti) are characterized by the highest microhardness. For 
these materials, the average scratch width at various forces was minimal. During tribological tests, the best frictional 
characteristics were recorded for stainless steel in combination with experimental alloy 2.22, and for the friction pair 
“titanium alloy VT3-1 — hard alloy 2.23”. The friction of this combination of materials was characterized by low 
friction coefficients with a low level of fluctuations, minimal wear of samples, and changes in the initial microrelief of 
their surfaces. 

Discussion and Conclusions. As a result of the research, the optimal friction pairs from the point of view of tribological 
interaction were established, specifically “titanium alloy VT3-1 — hard alloy 2.23” and “stainless steel 12X18N9T — 
hard alloy 2.22”. The frictional interaction for these combinations of materials is characterized by minimal volumetric 
wear, which will contribute to increasing the wear resistance of the tool in the areas of elastic contact on the front and 


rear surfaces. 
Keywords: hard alloys, wear resistance, stainless steel, titanium alloy. 
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Introduction. Various materials based on tungsten carbide, including hard alloys (HA), are widely used in 
many branches of modern production due to a number of advantages of their physical and mechanical properties [1-3]. 
The largest share (about 65%) of tool cutting materials (TCM) used in metalworking also belongs to HA that provide 
high cutting speeds when processing various structural materials [1, 3]. Currently, many directions related to improving 
the HA performance characteristics by various methods are being developed [4—6]. 

One of them is the creation of new binder compositions for carbide phases [7—9]. Experimental hard alloys 
(EHA) based on single-carbide alloy VK8 with various types of modified binders' have been developed at the Metal- 
Cutting Machines and Tools Department, Don State Technical University. Despite the fact that cobalt, because of its 
properties, is the most common binder for WC, the use of this metal is economically unprofitable due to its scarcity and 
high price. Experimental alloys are characterized by high values of the thermal entropy of the modified cobalt binder 
(Table 1) and, consequently, lower thermo-emf with respect to tungsten carbide, which increases the electrochemical 
stability of these materials [9-11]. Among the materials processed with a carbide tool, titanium-aluminum alloys and 
austenitic stainless steels can be distinguished. Due to a number of characteristics, these materials are used for the 
manufacture of parts in the most critical areas of mechanical engineering, including aerospace, nuclear, food, and 
medical production. Due to a number of characteristics, these materials are used for manufacturing parts in the most 
critical areas of mechanical engineering, including aerospace, nuclear, food and medical production. Cutting of titanium 
alloys and stainless steels is difficult due to their low thermal conductivity, high cutting forces, as well as unsatisfactory 
tribotechnical characteristics [12]. In this regard, the investigation of the features of the frictional interaction of newly 
developed TCM with the mentioned materials is an urgent task. Taking into account modern environmental and 
economic requirements for the organization of production, the processing of these structural materials takes place 
mainly in the mode of minimum output (MQL - minimum quantity lubrication) or without the use of lubricant-cooling 
process media (LCPM) [13, 14]. 

Then, the contact areas of the TCM and the workpiece material to be processed on the front and rear surfaces 
of the cutting tool can be considered as tribosystems operating in the friction mode without LCPM or in the mode of 
boundary friction. This does not exclude a significant proportion of metal contact. In this case, the operational 
parameters of the TCM, as an element of the friction couple, will be significantly effected by its tribotechnical 
characteristics in the dry friction mode. 

The presented paper is devoted to investigating tribotechnical characteristics of experimental hard alloys under 
dry friction on hard-to-cut materials, as well as to determining the optimal combination of “EHA — structural material” 
from the point of view of frictional interaction. This work is part of a complex of studies on physico-mechanical, 
tribological and cutting properties of experimental HA with modified cobalt binder. 

Materials and Methods. The following HA compositions were selected as objects of study of surface 
mechanical characteristics (Table 1). Tribological tests were carried out on square-section indenters (a=5 mm, Ra=0.1- 
0.12 um) of the three most promising compositions (2.21, 2.22, and 2.23) on tribometer T-11 (Poland) implementing a 
“finger-disk” friction scheme. During the experiments, a change in the friction force (F, N) was recorded depending on 
the friction path (ZL, m). Each experiment was repeated 3—5 times, the experimental results were processed using 
methods of dispersion analysis [0]. Titanium alloy VT3-1 and stainless steel 12X18H9T were selected as the material of 


the rotating disk (counterbody); the surface roughness of these samples was within R, 0.12—0.15 um. 


' Ryzhkin AA, Meskhi BCh, Bokov AI, et al. Hard alloy based on tungsten carbide (options). RF patent no. 2531332, 2014. (In Russ.) 


Machine building and machine science 


ep 
—_ 


http://vestnik-donstu.ru 


i) 


Advanced Engineering Research 2022. V. 22, no. 1. P. 50-56. ISSN 2687-1653 





Table 1 


EHA chemical composition and thermodynamic properties 


Entropy, J/mol-deg 
Carbide grade 


Composition 


we os 


po 
92.63% WC+7. _ Ld 20COr 
ow +5.03%Fet+0.82%Cu] 
92.38% WC+7.62% [3.6%Cot 
= +3.2%Fe+0.82%Cu] 


92.45% WC+7.55% [5.3%Co+l.43%Fe+0.82%Cu] 
91.95% WC+8.05% [5.65%Co+l.8%Mo+0.6%Ti] 
91.59% WC+8.41% [5.1%Co+2.7%Mo+0.61%Ti] 


2.24 90.62% WC+9.38% 
[3.34%Co+5.44%Mo+0.6%Ti] 
92% WC + [7.5-8]%Co, Fe<0.3% 
asic 


The studies were carried out at different sliding speeds and temperatures at constant load P=20 N. The mass of 





the samples was determined on LV 210-A balance. The roughness of the friction tracks on the counterbody after the 
experiments was measured on the Abris-PM7 profilometer (Russia). The dynamic microhardness of the EHA surfaces 
was determined using a scanning nanohardness tester NanoSCAN-O1 (Russia) through analyzing the thickness of 
scratches applied with different forces. The studies on worn surfaces of EHA samples were carried out on inverted 
ZEISS AxioVert. Al microscope. 

Research Results. The smallest width / of scratches applied by various forces F', belongs to alloys 2.22 and 


2.23, the binder in which was modified by Mo-Ti group (Table 2). 
Table Z 
Scratch width / under various forces F’, according to the results of EHA sclerometry 
h, wm 
Carbide grades 





Thus, these materials are characterized by the highest surface hardness at the micro-level, which implies better 
resistance to abrasive wear. It must be said that no noticeable changes in the surface microrelief were detected in these 
alloys under forces F,<15 N. The lowest microhardness according to the test results was demonstrated by composition 
2.24. 

To determine and compare the HA wear resistance, the mass loss of the indenter was measured for each value 


of friction path L, and then the volumetric wear of samples JV was determined (Fig.1, 2). 
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Fig. 1. Volumetric wear of EHA indenters under friction on titanium alloy VT3-1 at temperatures: a) 25° C; b) 300° C 


Under friction on a titanium alloy at different temperatures, the greatest volumetric wear JV was observed in 
alloy 2.22. The best wear resistance was demonstrated by composition 2.23 (Fig. 1). Under friction on stainless steel, 
the lowest values of parameter JV were recorded for composition 2.22 (Fig. 2). In this case, the greatest volume wear 
also belongs to the HA experimental compositions. At room temperature, alloy 2.21 demonstrates the highest wear 


intensity, and, when the friction zone is heated, the maximum wear values are fixed for composition 2.23. 


AV, mm AV, mm? 


1.0 


0.8 





200 400 600 800 1000 Lm 200 400 600 800 1000. L,m 
a) b) 


Fig. 2. Volumetric wear of EHA indenters under friction on stainless steel 12X18N9T at temperatures: a) 25° C; b) 300°C 


The surfaces of indenters made of alloys that have demonstrated maximum volume wear are characterized 
either by the predominance of worn areas, or the original surface of the material is preserved only in the form of 
individual rare fragments. Figure 3 shows pictures of the surfaces of indenters made of alloys 2.23 and VK8 after 


friction on 12X18N9T steel at a temperature of 300° C for friction path L=600 m. 


The surface of the more wear-resistant alloy 2.22 has a spotted structure with a predominance of initial unworn 
areas (Fig. 3 a). The surface of the base alloy VK8 is characterized by a large scale of destruction, rare fragments of the 
initial surface are discretely located, their total area is much smaller (Fig. 3 5). 

The frictional interaction of friction couples was estimated through comparing the average values of friction 


coefficient f,, and its fluctuations without taking into account the run-in stage I (Fig. 4 a, b). The comparison of the 
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friction coefficient fluctuations was carried out by standard deviation o,, of this parameter from /,, at the steady-state 


friction stage II. 





200 um 
a) b) 


Fig. 3. Comparison of surfaces of indenters of alloys 2.22 (a) and VK8 (db) after friction on 12X18N9T stainless steel 
(T=300° C, v=0.3 m/s): 1 — worn surface areas; 2 — fragments of the initial surface 


The surface of the more wear-resistant alloy 2.22 has a spotted structure with a predominance of initial unworn 
areas (Fig. 3 a). The surface of the base alloy VK8 is characterized by a large scale of destruction, rare fragments of the 
initial surface are discretely located, their total area is much smaller (Fig. 3 5). 

The frictional interaction of friction couples was estimated through comparing the average values of friction 
coefficient f,, and its fluctuations without taking into account the run-in stage I (Fig. 4 a, b). The comparison of the 
friction coefficient fluctuations was carried out by standard deviation o,, of this parameter from /,, at the steady-state 


friction stage II. 
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Fig. 4. Dependence of friction coefficients fon path L under friction on VT3-1 alloy (T=25° C, v=0.3 m/s): 
a) 2.23; b) VK8: I — run-in stage; II — stable friction stage 


In the case of friction on a titanium alloy at different temperatures, the average friction coefficients and their 
standard deviations for all EHA were higher than for the base grade VK8. However, the lowest values of these 
parameters among the experimental compositions belong to alloy 2.23 (Fig. 4 a). 

Under friction on steel 12Kh18N9T, both at room temperature and with heating, the largest values of 
parameters f,, and o,, belong to VK8 base alloy (at 25° C: f,,=0.72 and o,,=0.048; at 300° C: f,,=0.68 and o,,=0.032). In 


this series of experiments, the minimum coefficients of friction and its fluctuations were fixed for composition 2.22. 
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Under friction without heating for this material, the values of the estimated parameters were /,,=0.44 and o,,=0.025, and 
at 300° C, they increased to f,=0.57 and o,,= 0.029. 

The average roughness of the friction tracks on the counterbodies was measured and compared for the 
maximum value of the friction path. Under friction on steel 12Kh18N9T, the smallest values of this parameter belong to 
alloy 2.22. At a temperature of 25° C, the roughness was Ra» 7.=4.12 um, at 300° C — Ra) ».=5.12 um. The highest 
roughness values were recorded during friction of VK8 base alloy. At room temperature, the value of this parameter 
was Rapgs=5.07 um, with heating of the friction zone — Ragrs=5.95 um. High roughness values indicate that the 
frictional interaction of this material and stainless steel was accompanied by larger-scale destruction under the 
formation, and destruction of adhesive and cohesive seams on the surface of the counterbody material. 

In a series of experiments with titantum alloy VT3-1, the best indicators of counterbody roughness were 
recorded after friction of the experimental composition 2.23 (Raz 3=3.35 um at 25° C and Ra) 3=4.54 um at 300° C). 
The greatest surface damage was obtained for the samples after frictional interaction with alloy 2.22. The surface 
roughness during friction without heating was Ra» ».=6.88 um, with heating — Ra) ».=8.07 um. 

Discussion and Conclusions. As a result of the study of the tribological characteristics of experimental hard 
alloys under friction on hard-to-cut materials, the best combinations in terms of the frictional interaction of a pair of 
materials have been established. For stainless steel 12Kh18N9T, the best tribological parameters were recorded under 
friction in combination with alloy 2.22, for titantum alloy — with composition 2.23. The friction process for these 
combinations of materials at different temperatures is characterized by minimal volumetric wear, low coefficients of 
friction, and a smaller scale of destruction of the surfaces of both hard-alloy indenters and counterbodies made of 
structural materials. 
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Introduction. The significance of machine learning under the conditions of digital transformation of industry, and 
methods of implementing deep learning to provide the performance of trust management systems are considered. The 
necessity of using convolutional artificial neural networks for deep machine learning is determined. Various 
technologies and architectures for the implementation of artificial neural networks are briefly considered; a comparative 
analysis of their performance is carried out. The work objective is to study the need to develop new approaches to the 
architecture of computing machines for solving problems of deep machine learning in the trust management system 
implementation. 

Materials and Methods. In the context of digital transformation, the use of artificial intelligence reaches a new level. 
The technical implementation of artificial neural systems with deep machine learning is based on the use of one of three 
basic technologies: high performance computing (HPC) with parallel data processing, neuromorphic computing (NC), 
and quantum computing (QC). 

Results. Implementation models for deep machine learning, basic technologies and architecture of computing machines, 
as well as requirements for trust assurance in control systems using deep machine learning are analyzed. The problem of 
shortage of computation power for solving such problems is identified. None of the currently existing technologies can 
solve the full range of learning and impedance problems. The current level of technology does not provide information 
security and reliability of neural networks. The practical implementation of trust management systems with deep 
machine learning based on existing technologies for a significant part of the tasks does not provide a sufficient level of 
performance. 

Discussion and Conclusions. The study made it possible to identify the challenge of the computation power shortage 
for solving problems of deep machine learning. Through the analysis of the requirements for trust management systems, 
the external challenges of their implementation on the basis of existing technologies, and the need to develop new 


approaches to the computer architecture are determined. 
Keywords: deep machine learning, processor, trust system, information security, computer, artificial intelligence. 
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Introduction. Over the past 9 years, the fourth industrial revolution has been taking place on an increasingly 
large scale in the world, including in Russia. One of its key components 1s digital transformation affecting all aspects of 
economic life — from a large-scale industrial production to the service sector, science, education, and households. 
Under these conditions, the use of the machine in tasks with which a person does not cope or copes worse than the 
machine, significantly expands. If earlier it was about the mechanization of manual labor and production automation, 
now a human being can be replaced by a machine in solving the problems of data processing, analysis, forecasting and 
management of various systems: equipment, engineering processes, industrial enterprises, retail chains, etc. 

The practical basis for replacing a human being with a machine in certain areas of intellectual activity is the 
use of artificial intelligence. The creation of a “strong” artificial intelligence is a task of the future, associated with the 
need to create and develop new technologies, as well as solving significant ethical problems. Currently, machine 
learning systems of varying complexity are available for application, representing a step towards creating a “strong” 
artificial intelligence. 

The global market for machine learning systems is expanding rapidly. In 2020, its volume amounted to $11.33 
bln, in 2021, it grew to $15.50 bln, and by 2028, it will reach $152.24 bln, showing an average annual growth of 
38.6%. 

The scope of machine learning systems is very large and includes marketing and trade, banking, industrial 
production, medicine, etc. Machine learning systems are most in demand in the following industries: 

— robotics for the intellectualization of industrial and service robots, including collaborative ones; 

— automated control systems for processes and enterprises; 

— production process control systems; 

— supply chain and customer relationship management systems; 

— executive production systems; 

— production analytics systems for process equipment; 

— business intelligence systems, etc. 

Complex process systems, e.g., implemented in the designs of machine tools with real-time numerical control, 
currently cannot be equipped with control systems suitable for machine learning. This requires computing power 
capable of performing significant computations in tens of microseconds. Such a problem cannot be solved with the help 
of modern technical means. Therefore, in most cases, the process of machine learning of the control system is carried 
out on computers separated from it without time limits, and then the learning results are transferred to the control 
system in the form of recommendations, instructions on operating modes, tool changes, verification intervals, etc. 

Machine learning methods conditionally correspond to the types of inferences that underlie them: induction, 
deduction, and traduction. The case-based supervised learning method, when large amounts of data prelabeled by a 
human operator are loaded into the machine, corresponds to induction. The unsupervised learning method, when the 
machine itself must find patterns in the data, identify patterns, arrange and structure the data, corresponds to induction 
and traduction. The expert method based on the use of specified patterns and patterns for data processing correspond to 
deduction and traduction. Traduction is implemented mainly through the use of transfer learning, based on the 
application to a given task of knowledge gained in solving another task. 

Machine learning uses various technologies and mathematical models. The model of artificial neural networks 
(ANN), built by analogy with biological neural networks, 1.e., networks of nerve cells of a living organism, has the 
greatest potential for development. ANN is a system of interconnected and interacting artificial neurons implemented in 
the form of processors, processor elements in the form of accelerators or coprocessors under the control of a central 


' Machine Learning Market, 2021-2028. Hardware & Software IT Services Market Research Report. 2021. P. 160. URL: 
https://www.fortunebusinessinsights.com/infographics/machine-learning-market- 102226 (accessed: 06.11.2021) 
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processor. ANN neurons are located in levels (layers). The first level corresponds to receiving, processing input data, 
and passing them to the next level. Intermediate levels are hidden, their task is to process incoming data and transfer it 
to the last (output) level. A neural network may have several hidden levels interspersed with levels where logical, 
mathematical, and other transformations are performed. From level to level, the data is processed, at each subsequent 
level, the relationships of the previous one, are identified. Such a multi-level ANN has great potential and can be used 
to implement deep machine learning [1-3]. 

Deep machine learning is an ANN design method using multilayer filters to extract and model features from a 
set of input data*. Such learning can be supervised or unsupervised. It is also possible to use deep machine learning for 
expert systems. 

The technical implementation of artificial neural systems with deep machine learning is based on the 
application of one of three basic technologies: high-performance computing with parallel data processing, 
neuromorphic and quantum computing”. 

Materials and Methods 

High-performance computing. High-performance computing with parallel data processing is implemented 
through hybrid computing systems, 1.e., systems with a heterogeneous hardware computing structure, including a 
central processing unit (CPU) and an additional computing module in the form of an accelerator or coprocessor. 
Depending on the processors used for parallel data processing, hybrid computers have one of four architectures: 

1. Graphic processing unit (GPU) based architecture. The most common solutions are graphics accelerators 
that expand the computing capabilities of the central processing unit of a computer system. The latest advances in this 
field are NVIDIA Tesla V100 graphics accelerators, providing 120 TFLOPS performance for deep machine learning 
tasks, i.e., 1.21014 floating-point operations per second’*.This is 500-1000 times higher than the performance of an 
ordinary personal computer (PC). It should also be taken into account that the specified performance is provided when 
solving problems that require significant computing power and significant time costs, but not when working in real 
time. Currently, GPU-based architecture is the most accessible one. In particular, to implement a system with limited 
computing power, it is enough to have a video card with a NVidia graphics processor on a PC that implements the 
CUDA parallel computing hardware and software architecture. Along with CUDA, GPGPU technologies that use the 
graphics processor of a video card for computer graphics to perform mathematical calculations include AMD 
FireStream technology (for ATI graphics processors). The global GPU market is currently around $26 bln and is 
growing at a rate of up to 34% per year’. 

2. Architecture based on field-programmable gate arrays (FPGA) — semiconductor devices that can be 
reprogrammed and change the topology of connections in use. The rated performance of these devices is relatively low 
— about 20 TFLOPS, however, the efficiency of using computing power is the highest among all the considered 
architectures. It is 6-7 times higher than that of graphics accelerators. The high efficiency of FPGA is due to the 
flexibility and speed of adjustment to the computational tasks being solved. According to Grand View Research, in 
2020, the global FPGA market amounted to $ 9.85 bln, the expected market growth rate for the period up to 2027 is 
9.7% per year’. 


* Glek P. Deep Learning: short tutorial. neurohive.io. URL: https:// neurohive.io /ru/osnovy-data-science/glubokoe-obuchenie-deep-learning-kratkij- 
tutorial/ (accessed: 06.11.2021) (In Russ.) 

> Kak ~~ sokratit' izderzhki pri ispol'zovanii II. | HitachiVantaraCorporation. URL: _ https://hitachi.cnews.ru/articles/2021-06- 
14 kak_sokratit_izderzhki_pri_ispolzovanii (accessed: 07.11.2021) (In Russ) 

* Kak sokratit' izderzhki pri ispol'zovanii II. Op. cit. 

> Global Graphics Processing Unit (GPU) Market Insights and Forecast to 2027. QyYResearch. 2021. P. 116. URL: 
https://reports.valuates.com/market-reports/QY RE-Auto-25 V3358/global-graphics-processing-unit-gpu (accessed: 10.11.2021) 

° Field Programmable Gate Array Market, 2020 — 2027. Grand View Research, 2020. P. 130. URL: https://www.grandviewresearch.com/industry- 
analysis/fpga-market (accessed: 11.11.2021) 
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3. Architecture based on the advanced special integrated circuit (ASIC). Due to the narrow specialization of the 
computational problems to be solved, they can be much simpler, cheaper, and more compact. ASIC performance can 
reach 1000 TFLOPS, but the efficiency of using computing power, e.g., the number of recognized images, is about 2 
times lower than in graphics accelerators. The growth rate of the ASIC market is significantly lower than that of the 
GPU market. In 2020, according to Global Industry Analysts, the global ASIC market amounted to $17.3 billion, the 
expected average annual market growth until 2027 is 7.7%’. 

4. Architecture based on single-chip accelerators (SoC). “System on a chip” SoC 1s a fully functional electronic 
device that has a motherboard, processor, and other components required for operation, placed on a single integrated 
circuit. SoC are common in mobile computers (smartphones), single board computers, and other embedded systems. At 
the same time, SoC have a significant potential use as part of hybrid computers. In addition, solutions are possible in the 
form of a single-chip SoC assembly with FPGA elements (Xilinx® Versal architecture for adaptive computing). The 
SoC market is currently very large, and it was worth $79.7 bln in 2020. For the period up to 2027, the market is 
projected to grow by 4.4% per year, by 2027, the market volume will reach $107.4 bln $”. 

Figures 1 and 2 show data on the actual and projected growth of the global market for chips for deep machine 


learning, prepared by Omdia’® consulting company. 
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7’ ASIC - Global Market Trajectory & Analytics April 2021. Global Industry Analysts, Inc.; 2021. P. 185. URL: 
https://www.researchandmarkets.com/reports/5 14093 9/asic-global-market-trajectory-and-analytics (accessed: 11.11.2021) 

* Michael Feldman. Xilinx Unveils Its Most Ambitious Accelerator Platform. 2018. URL: https://www.top500.org/news/xilinx-unveils-its-most- 
ambitious-accelerator-platform/ (accessed: 10.11.2021) 

” System-On-A-Chip (SoC) - Global Market Trajectory & Analytics. Report, April 2021. Global Industry Analysts, Inc.; 2021. URL: 
https://www.researchandmarkets.com/reports/28323 1 6/system-on-a-chip-soc-global-market-trajectory (accessed: 11.11.2021) 

'° Joshi A. Deep Learning Chipsets Report — 2020. Omdia Marke, 2020. URL: https://omdia.tech.informa.com/products/deep-learning-chipsets- 
report---2020 (accessed: 11.11.2021) 

'! Joshi A. Deep Learning Chipsets Report — 2020. Op cit. (accessed: 12.11.2021) 
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The analysis shows that against the background of the general growth of the market, ASIC have the greatest 
prospects, GPU and SoC will also retain significant positions. For a fairly long term, GPU accelerators do not have an 
adequate replacement for solving complex problems, including in the learning process, and SoC are indispensable for 
mobile implementations of deep machine learning systems, as well as for parallel computing to offload the central 
processor. Rejection of FPGA in deep machine learning systems seems optional, although the share of such processors 
is likely to be relatively small. 

One of the most promising directions in the development of hybrid computers is the use of tensor and other 
specialized coprocessors, such as machine vision processors. Such coprocessors can be based on the most common and 
high-performance ASIC, as well as FPGA or GPU. The difference between coprocessors and accelerators is in the 
degree of integration with the central processor. The central processor translates control instructions to the accelerator 
through a special memory area. The coprocessor monitors the flow of machine code instructions from RAM to the CPU 
and intercepts instructions appropriate to its functional purpose, such as tensor transform tasks, pattern recognition, etc. 
To solve large-scale problems that require long-term distributed computing, it is advisable to use an accelerator; for 
frequent and repeated execution of simple computing tasks, use a coprocessor, in which the central processor is not 
loaded and does not slow down data processing. 

The operational properties of computers depend significantly on the architecture used. For solving deep 
machine learning problems, where up to 80% of computing power is spent, machines based on graphics accelerators are 
best suited. They have high performance in solving complex tasks that require a significant investment of time, have 
high flexibility and maximum calculation accuracy, but low relative performance. Specifically, for NVidia processors, it 
is 1.3-1.8 GOPS/W. ASIC-based computers have the highest absolute and relative performance. For neuIBM 
processors, e.g., the relative performance 1s 254 GOPS/W [4]. However, such machines have low flexibility and limited 
accuracy, So, it is advisable to use them when solving typical, e.g., matrix or tensor transformations, repetitive or multi- 
threaded tasks in real time. 

FPGA-based computers have high parameters of flexibility, accuracy, absolute and relative performance. For 
Tegra TX1 processors, e.g., the relative performance is 70 GOPS/W. However, such machines have a relatively high 
cost, so it is advisable to use them for scientific purposes, when only a few computers of a given configuration are 
required, as well as for developing the architecture of mass-produced ASIC and SoC processors. 

Despite the improvement in the architecture of computers, the potential for growth in capacities for high- 
performance computing will soon be exhausted. The number of transistors on a chip over the past 5 years has increased 
by about 12 times [5], and the amount of calculation in the process of machine learning has increased by 150 thousand 
times’. 

Neuromorphic computing. A possible way to address the lack of computing power for artificial neural 
systems with deep machine learning is to use neuromorphic computing and related chips. A neuromorphic chip is a 
processor based on the principles of the human brain. Such a chip simulates the work of neurons and their processes — 
axons and dendrites, which are responsible for the data transmission and perception. Connections between neurons are 
formed by synapses — special contacts through which electrical signals are transmitted. 

Some of the best-known developments in this area include IBM TrueNorth neuromorphic processors and Intel 
Loihi processors. They use an asynchronous cluster architecture and a convolutional neural network model — a 
unidirectional multilayer network with alternating convolutional and subsampling layers. TrueNorth processor is based 


on 28 nm technologies, Loihi — on 14 nm [6]. 


"Thompson NC, Greenewald K, Lee K, et al. The Computational Limits of Deep Learning. arXiv preprint arXiv:2007.05558. 2020 (accessed: 
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TrueNorth NS16e-4 multiprocessor system, containing 100 mln neurons and designed to work with networks 
for deep machine learning, was introduced by IBM in 2018 [7]. Each chip contains 1 mln digital neurons and 256 mln 
synapses enclosed in 4096 synapse nuclei; power consumption of each chip is 70 mW. 

Loihi processor, introduced in 2017, contains 131,000 artificial neurons and 131 mln synapses. In 2019-2020, 
Intel introduced two products based on Loihi — PohoikiBeach and PohoikiSprings. PohoikiBeach computing system, 
which includes 64 Loihi processors, has a total of 8.32 mln neurons and 8.32 bln synapses. PohoikiSprings Computing 
System Includes 768 Loihi processors, 100 mln neurons, and 100 bin synapses’. 

In Russia, work on the creation of neuromorphic processors has been going on for several years. In 2020, 
Motive Neuromorphic Technologies created the Altai neurochip [8]. The processor technology standard is 28 nm, 
power consumption is about 0.5 W, the crystal area is 64 mm/ (for comparison: TrueNorth - 430 mm’, Loihi - 60 mm’). 
It has 131 thousand neurons, between them, there are 67 mln connections. 

To assess the quality of neuromorphic processors, the following is used: 

1. Absolute performance indicator. This is the number of billions of synoptic operations performed per second 
(GSOPS). 

2. Energy efficiency indicator. This is the number of picojoules of energy expended in performing one synaptic 
operation (pJ/SOP). 

TrueNorth processor has a performance of 58 GSOPS and an energy efficiency of 26 pJ/SOP*. Similar power 
efficiency (23.7 pJ/SOP [9]) is provided by Loihi processor. 

The only competitor of neuromorphic processors in the implementation of neural networks with deep machine 
learning in the midterm (8—12 years) is hybrid computers with ASIC coprocessors. Such processors have a lower but 
comparable performance of synaptic operations and higher energy efficiency. In particular, ASIC processor described in 
[10] provides a synaptic performance of 8.7 GSOPS and an energy efficiency of 15.2 pJ/SOP. 

The global market for neuromorphic chips is relatively new and therefore small. In 2020, its volume amounted 
to only $22.5 mln. At the same time, the growth rate of the market is very high. By 2026, the market will grow to 
$333.6 mln, which corresponds to an average annual growth of 47.4%". 

Sometimes neuromorphic chips are understood as all types of processors that externally reproduce the work of 
neurons, regardless of the internal structure of the technical device, which may not correspond to the nature of the 
interaction of neurons. Such processors used to build artificial neural networks are properly called neural. Along with 
neuromorphic chips, neural processors also include processors (chips) with tensor and other specialized coprocessors 
for machine vision, speech recognition, etc. The world market for neural processors currently stands at $2.3 bln, and, by 
2027, it will grow to $10.4 bln, thus, the average growth will be 24.2% per year’®. 

Quantum computing. In the long term, the development of quantum computing can become a means of 
eliminating the shortage of computing power. Quantum computing solves problems through manipulating quantum 
objects: atoms, molecules, photons, electrons, and specially created macrostructures. Manipulations of quantum objects 
provide using: 

— quantum superposition, which manifests itself in the ability of quantum systems to simultaneously be in all 


possible states; 


'° Intel Scales Neuromorphic Research System to 100 Million Neurons. Intel, 2020. URL: https://newsroom.intel.com/news/intel-scales- 
neuromorphic-research-system-100-million-neurons/?utm_source=ixbtcom#gs.7oc6iw (accessed: 12.11.2021) 

'* Neurochip “Altai”. motivnt.ru. URL: https://motivnt.ru/neurochip-altai/ (accessed: 10.11.2021) (In Russ.). 

'S Neuromorphic Chip Market - Growth, Trends, COVID-19 Impact, and Forecasts (2021 - 2026). Mordor Intelligence, 2020. URL: 
https://www.mordorintelligence.com/industry-reports/neuromorphic-chip-market (accessed: 11.11.2021) 

'© Neuromorphic Chips - Global Market Trajectory & Analytics. Global Industry Analysts, Inc. 2021. 118 p. URL: 
https://www.researchandmarkets.com/reports/4805280/neuromorphic-chips-global-market-trajectory-and (accessed: 10.11.2021) 
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— quantum entanglement, which manifests itself in a strong relationship between the parameters of specially 
prepared quantum systems. 

Devices for quantum computing are usually divided into two large classes [11]: general purpose quantum 
computers and quantum simulators. The former, like central processing units, can solve any algorithmic problem, and 
quantum simulators are analog computers for solving highly specialized problems. 

Technologies for creating universal quantum computers are currently at the stage of formation. The created 
computers demonstrate “quantum superiority” in solving certain problems, but so far, they cannot be used to form 
artificial neural networks with deep machine learning. Companies that are most active in creating a quantum computer 
include: 

1. Google. In 2018, a 72-qubit Bristlecone quantum processor was built, in 2019, a more accurate 53-qubit 
Sycamore quantum processor was built. 

2. Intel. A 49-qubit TangleLake superconducting quantum chip was built in 2018. 

3. IBM. In 2017, a 50-qubit quantum processor was created and tested, in 2019 — the world's first 
commercial 20-qubit quantum computer IBM Q SystemOne, etc. 

The only adiabatic quantum computer on the market 1s D-WaveSystems, available in 16 to 2000 qubits, 
arranged in clusters of 8 qubits each. 

The field of quantum simulators is also rapidly developing. One of the most complex simulators of this type is 
a 2017 joint development of the University of Maryland and the National Institute of Standards and Technology (USA). 
This 53-qubit simulator uses cold ytterbium ions as qubits. A similarly capable 51-qubit quantum simulator based on 
rubidium atoms was developed by a group of scientists at Harvard University and Massachusetts Institute of 
Technology. 

A number of projects developing quantum computing technologies are also being implemented in Russia. In 
particular, for several years now, the development of a superconducting processor has been carried out by scientists 
from a consortium, which includes National University of Science and Technology (MISIS), Osipyan Institute of Solid 
State Physics of the Russian Academy of Sciences (ISSP RAS), Institute of Solid State Physics, Bauman Moscow State 
Technical University (MBSTU), Dukhov Automatics Research Institute (VNIIA), and other organizations. To date, the 
consortium has debugged the technology for manufacturing superconducting two-qubit circuits, experimentally 
characterized and demonstrated two-qubit logic gates that perform quantum entanglement, which is required for the 
operation of quantum processors. The reliability of logical operations is in the range of 85—95%. 

In 2020, $675 mln was invested in quantum computing in the world, which is more than 3 times the investment 
amount in 2019 ($211 mln). In 2021, the volume of investments in quantum computing exceeds $800 mln [12]. 

Trust management systems. One of the basic requirements for management systems, including production 
ones, is to provide them with the required level of trust. According to GOST R 54583-2011 “Information technology. 
Security techniques. A framework for IT security assurance. Part 3. Analysis of assurance methods”, the purpose of 
providing credibility is to create confidence in the reliable functioning of the product under given conditions. To 
provide this, the information system must have the following operational properties [13]: 

— functional reliability, 1.e., the ability to perform its function with a given reliability, which in turn is 
normalized by the number of failures, the error and repeatability of the calculation results; 

— information security, 1.e., the ability to provide a given level of confidentiality, availability and integrity of 


information: stored, transmitted, received, and processed during the operation of the system. 
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The subject of this research is control systems that relate to information. Therefore, the above requirements for 
performance properties are also valid for them. However, management systems have their own specifics in the 
definition of trust. A trust control system must have: 

— the ability to control, e.g., a robot, a machine tool, an enterprise, etc., according to a given number of 
parameters with a specified reliability, which is regulated by the number of failures, error, and repeatability, and with a 
given performance, which in turn is regulated by the time of data processing and execution of control commands; 

— the ability to control the elements, structure and processes of the system at the hardware and software 
levels to provide information security. 

If the control system has the function of deep machine learning, then the fulfillment of the first of these 
requirements imposes severe restrictions on the means of technical implementation used. This should be the optimal 
computer for the formation of a convolutional neural network with high parameters of performance, accuracy, and 
calculation reliability. 

If we do not consider the option of using quantum processors, the full-featured implementations of which are 
not yet available, then neither hybrid computers based on all the considered architectures, nor computers based on 
neuromorphic processors fully comply with the first requirement. Computers based on ASIC and neuromorphic 
processors do not provide high accuracy and reliability, and hybrid computers with GPU or SoC accelerators are not 
optimal for real-time operation, including impedance. 

A certain compromise is provided when using FPGA-based hybrid computers, however, such machines have a 
high cost in mass production, significantly lower performance than ASIC machines, and significantly less complex 
computing capabilities than GPU machines. Another compromise option is the simultaneous use of a CPU with a 
graphics accelerator to solve complex tasks in the process of machine learning and tensor or other highly specialized 
ASIC-based coprocessors for real-time data processing. 

The second requirement, although technical in content, in practice acts as an economic one. The 
implementation of a control system for process equipment with deep machine learning is possible only through a 
convolutional neural network, the control of which from the outside 1s not possible. Information security can be only 
provided that the main part of the computer will be created by domestic manufacturers who have been certified in the 
field of information security. 

At present, the main part of control systems in Russia is built on the basis of foreign microelectronic 
components. The share of such components exceeds 85% [14]. Providing information security in the case of using 
imported components in computers does not have an unambiguous solution and depends on the structure of the created 
artificial neural network and the order of its use. In particular, when using hybrid computers, information security can 
be significantly improved through localizing data transfer between the central processor and the accelerator or 
coprocessor. 

Research Results. The analysis of deep machine learning models, basic technologies, and architecture of 
computers, as well as the requirements for providing confidence in control systems using deep machine learning, allows 
us to draw the following conclusions: 

1. There is an objective problem of lack of computing power for solving problems of deep machine learning. 
None of the currently existing technologies can solve the full range of training and impedance problems. 

2. Since deep machine learning is implemented on the basis of a model of convolutional neural networks, their 
external control to provide information security and reliability of work is not possible. The only option is developer 
control, which also has limited capabilities. This determines the need for the production of processors required for ANN 


in Russia. 
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3. The practical implementation of trust control systems with deep machine learning based on existing 
technologies for a significant part of the tasks in real time cannot be provided, for the other part of the tasks, such an 
implementation is associated with a significant drop in performance. 

4. The performance gain of trust control systems can be based on improving the architecture of hybrid 
computers, including the simultaneous use of processors of different architectures that are optimal for solving the 
corresponding problems of analysis and control. 

Discussion and Conclusions. This paper analyzes and discusses the relevance and implementation of machine 
learning in the context of digital transformation of industry. The scientific problem covered in the work is in the 
insufficient development of the technical level of modern computers to provide high performance of algorithms based 
on deep machine learning. Attention is drawn to the problem of information security, which is one of the prerequisites 
for the development of domestic processors for ANN. Based on the analysis of the requirements for trust control 
systems, the objective difficulties of their implementation based on existing technologies and the need to develop new 


approaches to the architecture of computers are determined. 


References 

1. Zelensky A, Semenishchev E, Alepko A, et al. Using neuro-accelerators on FPGAs in collaborative robotics 
tasks. SPIE Optical Instrument Science, Technology, and Applications II. 2021;11876:1187600. 
https://doi.org/10.1117/12.2600582 

2. Zelenskii AA, Pismenskova MM, Voronin VV. Control of Collaborative Robot Systems and Flexible 
Production Cells on the Basis of Deep Learning. Russian Engineering Research. 2019;39:1065—1068. 
http://dx.doi.org/10.3103/S1068798X 19120256 

3. Voronin VV, Sizyakin RA, Zhdanova M, et al. Automated visual inspection of fabric image using deep 
learning approach for defect detection. Automated Visual Inspection and Machine Vision IV. 2021;11787:117870 
http://dx.doi.org/10.1117/12.2592872 

4. Phi-Hung Pham, Jelaca D, Farabet C, et al. NeuFlow: Dataflow Vision Processing System-on-a-Chip. In: 
Proc. IEEE 55th International Midwest Symposium on Circuits and Systems (MWSCAS). 2012. P. 1044-1047. 


http://dx.do1.org/10.1109/MWSCAS.2012.6292202 
5. Shuremov EL. Whether it is worth being fond of big data? Accounting. Analysis. Auditing. 2020;7:17—29. 


https://do1.org/10.26794/2408-9303-2020-7-2-17-29 

6. Jing Pei, Lei Deng, Sen Song, et al. Towards artificial general intelligence with hybrid Tianjic chip 
architecture. Nature. 2019;572:106—111. http://dx.doi.org/10.1038/s41586-019-1424-8 

7. Modha D. TrueNorth: from zero to 64 million neurons. Open Systems. DBMS. 2019;3:8. 

8. Akopyan A, Sawada J, Cassidy A, et al. TrueNorth: design and tool flow of a 65 mw 1 million neuron 
programmable neurosynaptic chip. IEEE transactions on computer-aided design of integrated circuits and systems. 
2015;34:1537-1557. http://dx.doi.org/10.1109/TCAD.2015.2474396 

9. Mike Davies, Narayan Srinivasa, Tsung-Han Lin, et al. Loihi: A neuromorphic manycore processor with on- 
chip learning. IEEE Micro. 2018;38:82—99. http://dx.do1.org/10.1109/MM.2018.112130359 

10. Kim J, Koo J, Kim T, et al. Efficient synapse memory structure for reconfigurable digital neuromorphic 
hardware. Frontiers in neuroscience. 2018;12:829. http://dx.do1.org/10.3389/fnins.2018.00829 

11. Fedorov A. Kvantovye vychisleniya: ot nauki k prilozheniyam. Open Systems. DBMS. 2019;3:14. (In 
Russ.) 

12. Bobier J-F, Langione M, Tao E, et al. What Happens When ‘If Turns to ‘When’ in Quantum Computing? 
BCG Digital Transformation. 2021. 20 p. 


Information technology, computer science, and management 


ON 
nn 


http://vestnik-donstu.ru 


Advanced Engineering Research 2022. V. 22, no. 1. P. 57-66. ISSN 2687-1653 





13. Sabanov AG. Doverennye sistemy kak sredstvo protivodeistviya kiberugrozam. ZaSita informacii. Inside. 


2015;63:17—21. 
14. Kalyaev IA, Melnik EV. Trusted control systems. Mechatronics, Automation, Control. 2021;22:227—236. 


https://doi.org/10.17587/mau.22.227-236 (In Russ.) 


Received 27.12.2021 
Revised 17.01.2022 
Accepted 18.01.2022 


About the Authors: 

Zelensky, Alexander A., Director of the Institute of Digital Intelligent Systems, Moscow State University of 
Technology “STANKIN” (3a, Vadkovsky Lane, Moscow, 127005, RF), Cand.Sci. (Eng.), associate professor, 
ResearcherID, ScopusID, ORCID, Zelenskyaa@gmail.com 

Abdullin, Tagir H., lecturer of the Industrial Electronics and Intelligent Digital Systems Department, Moscow 
State University of Technology “STANKIN” (3a, Vadkovsky Lane, Moscow, 127005, RF), senior engineer, ScopusID, 
ORCID, everestultimate@yandex.ru 

Zhdanova, Marina M., Junior Research Scholar, Moscow State University of Technology “STANKIN” 
(3a, Vadkovsky Lane, Moscow, 127005, RF), ResearcherID, ScopusID, ORCID, mpismenskova@mail.ru 

Voronin, Viacheslav V., Associate Director of the Center for Cognitive Technologies and Machine Vision, 
Moscow State University of Technology “STANKIN” (3a, Vadkovsky Lane, Moscow, 127005, RF), Cand.Sci. (Eng.), 
associate professor, ResearcherID, ScopusID, ORCID, voronin_sl@mail.ru 

Gribkov, Andrey A., Director of the Analytical Center, Moscow State University of Technology “STANKIN” 
(3a, Vadkovsky Lane, Moscow, 127005, RF), Dr.Sci. (Eng.), professor, ResearcherID, ScopusID, ORCID, 


andarmo@yandex.ru 

Claimed contributorship 

A. A. Zelensky: basic concept formulation; research objectives and tasks setting; text preparation; formulation 
of conclusions. T.H.Abdullin and M.M. Zhdanova: conducting research; analysis of existing approaches. 


V.V. Voronin and A. A. Gribkov: analysis of the research results; the text revision; correction of the conclusions. 


All authors have read and approved the final manuscript. 


Advanced Engineering Research 2022. V. 22, no. 1. P. 67—75. ISSN 2687-1653 





INFORMATION TECHNOLOGY, COMPUTER 
SCIENCE, AND MANAGEMENT 





® Check for updates 


UDC 004.11.5 Original article 
https://doi.org/10.23947/2687-1653-2022-22-1-67-75 lie, o 





Comparison of machine learning models for coronavirus prediction 


Brou Kouame Amos ><, Ivan Smirnov , Mabouh Moise Hermann ol 
Peoples' Friendship University of Russia (RUDN) (Moscow, Russian Federation) 


bX] broukouameamos9@gmail.com 

Introduction. Coronavirus, also known as COVID-19, was first detected in Wuhan, China, in December 2019. It is a 
family of viruses ranging from the common cold to severe acute respiratory syndrome (SARS). The symptoms of such a 
virus are similar to those of a cold or seasonal allergies. Like other respiratory viruses, it 1s mainly transmitted through 
airborne droplets when coughing or sneezing. Therefore, the recognition of COVID-19 requires careful laboratory 
analysis, and the reduction of recognition resources is a major challenge. On 11 March, 2020, the World Health 
Organization (WHO) declared COVID-19, caused by SARS-CoV-2, a pandemic, as there had been an exponential 
increase in cases worldwide, and demand for intensive beds and related structures had far exceeded existing capacity. 
The first examples of this are the regions of Italy. Brazil registered the first case of SARS-CoV-2 on 02/26/2020. 
Transmission of the virus in this country shifted very quickly from imported cases to local and, finally, community 
missions, with the Brazilian federal government announcing national community transmission on 03/20/2020. As of 
March 23, in the state of SAo Paulo with a population of about 12 million people, where the Israelita Albert Einstein 
Hospital is located, 477 cases of the disease and 30 related deaths were registered, and on March 27, there were already 
1223 cases of COVID-19 with 68 concomitant deaths. To slow the spread of the virus in the state of Sao Paulo, 
quarantines and social distancing measures were introduced. One of the motivations for this challenge is the fact that, in 
the context of an extensive healthcare system with the possible limitation of SARS-CoV-2 testing, it is not practical to 
test every case, and test results can only be used in testing the target subpopulation. The study objective is to build a 
model based on machine learning that can predict the detection of SARS-CoV-2 from medical data. For this, various 
classification models of machine learning are compared, and the best one to predict coronaviruses is determined. The 
comparison is based on individuals in class 1, 1.e., those with a positive test. Therefore, it is required to determine the 
machine learning model with the best response and F1 score for class 1. 

Materials and Methods. An open-source data set from the Israelita Albert Einstein Hospital in Sao Paulo, Brazil, was 
taken as a basis. The following machine learning models were used for the study: RandomForests (RF), K-Nearest 
Neighbor (KNN), Support Vector Machine (SVM), Logistic Regression (LR), Decision Tree (DT) and AdaBoost (AB), 
as well as the 10-time cross-validation technique. Some machine learning performance measures, such as accuracy, 
recall, and F1 score were evaluated. 

Results. Out of a total of 5,644 people tested during the COVID-19 pandemic, 5,086 people tested negative and 
558 people tested positive. At the same time, support for machine vectors showed the best results in detecting 
coronavirus with a recall of 75 % and an F1 score of 60 % compared to models: Random drill, KNN, LR, AB, and DT. 
Discussion and Conclusions. It was found that when using AB algorithms, greater accuracy is achieved, but the 
stability of the LSVM algorithm is higher. Therefore, it can be recommended as a useful tool for detecting COVID-19. 


Keywords: COVID-19 detection, classification, machine learning models. 
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The coronavirus is a very severe acute respiratory syndrome caused by the SARS-COV-2 virus. This virus, 


1. Introduction 


which can infect humans or animals, was discovered in the Chinese region of Wuhan, more precisely in the province of 
Hubei, during the pneumonia epidemic of January 2020 [1,2]. It is therefore the seventh human coronavirus. To 
everyone's surprise, this virus spread worldwide, causing 318,599 deaths and 4,806,299 infected persons [3]. 

SARS-CoV-2, SARS-CoV and MERS-COV (Middle East Respiratory Syndrome Coronavirus) cause severe 
pneumonia with a mortality rate of 2.9 %, 9.6 % and 36 % respectively [4—6]. 

The other four viruses, namely OC43, NL63, HKU1, and 229E, are responsible for illnesses related to mild 
symptoms |7]. 

It should be noted that since the Covid-19 epidemic, there has been much speculation about the origin of this 
virus [8]. Some said that it was the result of work done in a laboratory. However, after studies conducted on genetic 
data, this hypothesis was dismissed [9]. Analysis and comparison with the genomes of previously known coronaviruses 
clearly show that SARS-COV-2 is different from other coronaviruses [8, 11]. The virus responsible for the coronavirus 
(SARS-COV-2) is similar to the SARS virus of bats [2]. Thus, the Covid-19 virus is believed to have originated from a 
bat coronavirus that became infectious to humans while acquiring genes specific to pangolin coronaviruses. It should be 
noted that the actual causes of Covid-19 are still unclear. 

The symptoms of Covid-19 are similar to those of seasonal flu. The disease is more severe in the elderly and 
in people who are vulnerable to certain chronic diseases. Patients with Covid-19 can have symptoms ranging from mild 
to severe. The most common symptoms are fever (83 %), cough (82 %) and breathlessness (31 %) [12]. In patients with 
pneumonia, the X-ray of the lungs shows numerous mottles and ground glass opacity [12, 13]. 

Gastrointestinal symptoms associated with patients with Covid-19 include vomiting, diarrhoea, and abdominal 
pain [12, 14]. 

We also see a decrease in lymphocytes and eosinophils, lower haemoglobin levels, and an increase in white 
blood cells and neutrophils [15-18]. 

The manifestation of Covid-19 in children is different from that in adults. In children, the symptoms are mild. 
However, in some children, we have seen severe and fatal cases [19-27]. 

Like all other viruses, Covid-19 is transmitted mainly by the respiratory route. Among these routes of 
transmission, we have droplet transmission, which is the most widespread [28, 29]. Other transmission routes exist, 
namely the faecal route, via saliva. Indeed, SARS-CoV-2 RNA was found in the stool of a patient with Covid-19 [31]. 
SARS-CoV-2 RNA can be detected on inanimate surfaces (door handles). People who have been in contact with these 
surfaces could be contaminated [29]. 

This model will make it possible to identify positive and negative cases from the dataset studied and the 
elements responsible for COVID-19. The proposed prediction model ensures that it tracks the results regarding this 
epidemic situation so that the huge economic losses, the spread of the community, the amount of detachment social 
gens can be detected and a precise decision can also be made accordingly. This method will allow government 


authorities to put in place preventive measures based on our future work to predict the onset of this disease in the future. 
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2. Data Resources and Methods 

The dataset used was uploaded to Kaggle. It is open source and available on this link 
kaggle.com/einsteindata4u/covidl19. This dataset contains anonymized data in accordance with best international 
practices and patient recommendations at the Israelita Albert Einstein Hospital in Sao Paulo, Brazil. This section 
describes the proposed approach and a detailed overview of the tasks. These tasks can help to understand and extract 
knowledge from COVID 19 data, which can help countries contain the spread of the virus, raise awareness, launch 
initiatives, determine if mitigation has a positive effect or not, identify other factors affecting the virus, etc. This will 
allow countries to prepare for what may happen in the near future. This could help save lives and alleviate the agony. 
Epidemiological information includes various characteristics of the case studied, including case identification, age, sex, 
target value, lymphocytes, leukocytes, monocytes, hco3, etc. 

2.1. Data Pre-processing 

In data analysis, the most important step is pre-processing. However, it is not clear what methods of pre- 
treatment the author used. This part must be completed. 

2.2. Data Transformation 

The data is transformed to be processed and stored in. xls for further processing. All data were normalized to 
have a mean of zero and a unit standard deviation. With a dataset containing 111 characteristics, data mining eliminated 
missing values (78 characteristics) and retained important characteristics (33). This exploratory analysis of the data also 
allowed us to identify two categories of characteristics, namely virus-related characteristics and blood-related 
characteristics. The target value is divided into two categories which are negative cases coded by 0 and positive cases 
coded by 1. 

The dataset from the Israelita Albert Einstein Hospital in Sao Paulo is divided into training and test data. 70 % 
of the data is used for predictive model training, and the remaining 30 % is used for testing. The objective of model 
training 1s to adapt the model using data from the training set. After the model is formed, the prediction models sound 
tested to evaluate performance in the test datasets. 

2.3. The Proposed Models 

This section describes the different machine learning models used in this paper. These models are: Random 
Drills (RF), K-plus Close Neighbors (KNN), Linear Support Vector Machine (SVM), Logistic Regression (LR), 
Decision Tree (DT), and AdaBoost (AB). 

Random Forest (RF) 

Random forests (RF) or random decision forests were first proposed in 1995. This is a general classification 
training method that tends to work better than traditional decision tree classification methods (Gangaie et al., 2019). 
Decision trees are the fundamental RF classifiers that vote for each of the forecasts, and the survival prediction is based 
on the majority voting method in each tree (Breiman, 2001). The accuracy of each tree and the independence of the 
trees from each other provide the reliability of the classification. We used 100 trees to predict two target classes, 
survival or death of patients with hepatitis. 

Nearest Neighbor (KNN) 

The K-Nest Neighbor (KNN) classifier is one of the most commonly used classification algorithms. This 
algorithm can be used in several applications. It saves all valid attributes and classifies new attributes according to their 
similarity dimension. KNN 1s a statistical recognition model method for detecting the different classes of a model. A 
tree data structure is used to determine the distance between the point of interest and the points in the training dataset. 
The attribute is classified by its neighbors. In the classification method, the value of k is always a positive integer 


closest to the neighbor. The nearest visions are selected from a set of classes or property values of the object. 
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Support Vector Machine (SVM) 

SVM-controlled learning method is used for classification and regression [29]. This algorithm is a relatively 
new approach and has performed well in recent years. The SVM classifier is based on linear classifiers and in the data 
separated by a row, the SVM isolates the objects in the specified classes. It can also identify and classify instances that 
are not supported by the data. The only extension of this algorithm is to perform a regression analysis to obtain a linear 
function, and another extension teaches to classify the elements to obtain a classification of individual elements. 

Logistic Regression Model (LR) 

Logistic regression is the corresponding regression analysis that should be performed when the dependent 
variable is dichotomous (binary). Like all regression analyses, logistic regression is predictive analysis. It is used to 
describe the data and explain the relationship between a dependent binary variable and one or more nominal, ordinal, 
interval or ordinal independent variables, report [30,31]. This approach assumes that the binary result follows a 
binomial distribution. 

Decision Tree (DT) Model 

The Decision Tree is a controlled learning method that is used to solve classification and regression problems, 
but it is more used to solve classification. This is a powerful classification method for disease prediction. This is a tree 
model where the internal nodes represent the characteristics of a data set, the branches represent the decision rules, and 
each leaf node represents a result. The decision tree consists of two nodes, a decision node and a leaf node. Decision 
nodes have multiple branches and are used to make a decision, while leaf nodes are the result of those decisions. 

Model AdaBoost (AB) 

AdaBoost, short for “Adaptive Boosting”, is the first boost algorithm proposed by Freund and Schapire in 
1996. Its goal is to turn weak predictors into strong predictors to solve classification problems. For classification, the 
final equation can be put under the heading below: 

F(x)= sign( Yin=1 9m fm@ ) (1) 
Where f,,, denotes the weak classifier m and m denotes the corresponding weight. AdaBoost can be used for face 
recognition, as it is a standard algorithm for detecting faces in images. AdaBoost is fast, requires no setup, and is simple 
and easy to program. Plus, it has the flexibility to be able to be combined with any machine learning algorithm. 

2.4. Evaluation of Performance Measures 

For the comparison of the different classification algorithms used in this paper, some metrics were evaluated. 
These are accuracy, recall, and Fl-score. These metrics are calculated based on true positives (TP), true negatives (TN), 
false positives (FP), and false negatives (FN). The standardized confusion matrix illustrates the relationship between 
classification results and predicted classes. The level of the classification performance is calculated by the number of 
samples correctly and incorrectly classified in each class. 


The accuracy is calculated based on the total number of correct predictions, defined as follows: 


TP+IN 
TP+FN+ITN+FP 


(2) 


Recall, or sensitivity, 1s the proportion of true positive predictions that have been correctly identified, defined as 


Accuracy = 





follows: 
Recall = — (3) 
TP+FN 
The F1 score is the harmonic mean of accuracy and recall, and it is calculated by: 
Score Fl = — (4) 
TP+5(FP+FN) 
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3. Result 


The objective of this paper is to compare the different models of machine learning for the detection of 
coronavirus. Our task was to find out which machine learning model has the best recall and fl-score for Class 1. The 
learning machine models used are: Radom drill, k-nearest neighbor, logistic regression, support vector machine, 
AdaBoost, and decision tree. Out of a total of 5,644 people tested for COVID-19, 5,086 people tested negative and 558 
people tested positive. The results of our study are presented in Figure 1 and Figure 3. These results show that the 
vector-machine gave better results with a recall of 75 % and an F1 score of 60 %. The different learning curves were 
also traced in order to understand the phenomenon of over-fitting and under-fitting Figure 2. Indeed, the learning curve 
is very well known to data scientists, the learning curve shows the efficiency and quality of learning of our machine 
learning model. Learning curves are widely used as a diagnostic tool in machine learning for algorithms that 
incrementally learn a training data set. This means that we increase our dataset by a certain step, and then we see the 


performance of our model. The model can be evaluated on the training dataset and on the exception validation dataset 


after each update during training, and it traces the measured performance. This can be represented as a curve. 
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Fig. 1. Classification report of different machine learning models 
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Fig. 3. Results of predictions from various machine learning techniques 
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Figure 3 shows the performance of the different machine learning algorithms according to the performance 
measures used in this paper. We see that for recall and Fl-score, LSVM outperforms the other machine learning models 
used, namely LR, KNN, RF, AB, and DT. For accuracy, LR is much better than the others. As for accuracy, we find 
that LR and AB performed better than the other models. In this paper, we chose recall and F1 score to measure the 
performance of the model. Recall allowed us to correctly identify the Covid-19 positive test subjects among all the real 
positive cases. As for the Fl score, we used it because we had an imbalance between different classes, 1.e., positive and 
negative cases. 

4. Discussion and Conclusion 

The data used in this paper was collected at the Israelita Albert Einstein Hospital in Sao Paulo, Brazil. After an 
exploratory analysis, two categories of characteristics were identified. These are the characteristics related to the virus and 
the characteristics related to the blood. Out of a total of 5,644 people tested with COVID-19, 5,086 people tested negative 
and 558 people tested positive. The results of this study clearly illustrated that in relation to our goal, machine vector 
support showed better results in coronavirus detection with a recall of 75 % and an F1 score of 60 %. This co-calculation 
was done with the other machine learning models, namely the Radom drill, the k-nearest neighbor, the logistic regression, 
the AdaBoost, and the decision tree. As such, this model can be useful for the diagnosis of COVID-19. However, it is 
possible to optimize the parameters of this model in order to improve its performance. 

After the analysis of the learning curve in Figure 2, we find that apart from the supporting sensor, other 
machine learning models can be studied for the detection of COVID-19. These include AdaBoost and k-nearest 
neighbor. Indeed, we find that if we perform a little more advanced optimization of the parameters of these models, they 
could be candidates for the diagnosis of COVID-19 because the difference between the learning score curve and the 


validation score curve would have reduced the model's ability to generalize. 
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